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SUMMARY 


Procedures  for  evaluating  the  probability  of  failure  of  structures  sub¬ 
jected  to  blast  and  shock  loadings  have  been  established.  For  simple  systems 
that  can  be  represented  by  a  slngle-degree-of-freedom,  a  procedure  Is  used 
that  computes  the  transition  probabilities  between  states  of  the  system  using 
conditional  probability  over  calculational  time  intervals.  The  system  and 
loading  parameters  that  are  random  must  be  adequately  characterized  and  the 
response  space  must  be  discretized  In  a  reasonable  way  to  represent  various 
states  of  response.  The  result  is  a  representation  of  the  distribution  of  the 
response  within  the  discretized  response  space.  All  other  probaDility  outside 
tne  discretized  response  space  is  lost  so  that  the  first  passage  probability 
outside  the  selected  bounds  of  response  can  be  obtained.  Various  schemes  to 
enhance  this  procedure  computationally  have  been  demonstrated.  The  procedure 
is  not  limited  to  elastic  systems  or  to  stationary  Input,  and  measures  of  the 
response  other  than  displacement  may  be  considered.  The  first  passage  of 
dissipated  energy  in  an  elasto-plastlc  system  was  specifically  addressed  and 
the  probability  of  failure  given  a  random  failure  level  has  been  computed. 

for  more  complex  systems  that  require  multiple-degree-of-freedom  repre¬ 
sentation,  a  proceaure  is  presented  that  expands  the  peak  response  in  a  Taylor 
series  in  terms  of  the  loading  and  system  parameters.  The  peak  response  is 
determined  by  conventional  deterministic  methods.  The  mean  and  variance  of 
the  response  at  each  location  are  determined  by  nunerically  approximating 
partial  difference  terms.  Given  the  mean  and  variance  of  the  failure 
criteria,  a  margin  of  survival  can  be  computed.  When  the  form  of  the 


distribution  of  the  margin  of  survival  Is  specified,  then  the  probability  of 
survival  and  Its  complement,  and  the  probability  of  failure  can  be  nunerlcally 
determined  by  integration  of  a  multivariate  distribution  function. 

An  additional  approach  for  complex  systems  has  been  demonstrated  that 
uses  the  relationships  between  parameters  in  a  numerical  sequence.  The  Markov 
property  allows  the  random  character  of  the  response  measures  and  their 
correlations  to  be  computed  at  each  time  step.  The  result  Is  the  evolution  of 
a  stochastic  process  with  time. 
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Section  II,  “Stochastic  Analysis  of  Mechanical  Systems,"  was  funded  under 
an  Air  Force  Office  of  Scientific  Research  (AFOSk)  grant,  program  element 
61101F,  and  has  oeen  published  elsewhere  as  part  of  AFWL-TR-82-123.  Sec¬ 
tions  III  and  IV  were  funded  under  program  elements  62601F  and  647IIF  respec¬ 
tively.  This  technology  is  being  transitioned  to  field  problems  concerned 
witn  random  loads  and  random  material  properties. 
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I.  INTRODUCTION 


1.  MOTIVATION  AND  OBJECTIVES 

Structures  are  designed  for  a  variety  of  purposes.  Some  function  as 
weather  shelters  (conventional  buildings)  while  others  protect  their  occupants 
from  more  severe  outside  environments  like  blast  and  shock  (protective  struc¬ 
tures).  Some  protect  tne  outside  environment  from  their  contents  (containment 
structures).  The  purposes  of  a  structure  are  usually  translated  into  perfor¬ 
mance  criteria.  The  goal  of  the  structural  engineer  is  to  design,  analyze, 
and  construct  structures  tnat  satisfy  the  performance  criteria.  Since  the 
inputs,  the  material,  and  tne  structural  performance  are  not  strictly  determi¬ 
nistic,  tne  engineer  must  oe  concerned  with  the  proDaoility  tnat  the  structure 
will  respond  to  all  applied  inputs  within  the  limits  set  by  the  criteria. 

This  proDability  is  called  reliability.  When  a  structure's  response  is 
designed  to  be  far  into  the  inelastic  regime,  as  in  the  case  of  most  protec¬ 
tive  structures,  the  consideration  of  reliability  becomes  botn  more  important 
and  more  difficult.  The  difficulty  arises  from  the  requirement  to  predict  the 
probabilistic  nature  of  trie  nonlinear  response. 

The  reliability  of  a  structure  snould  be  determined  in  any  comprehensive 
analysis  since,  in  practice,  most  inputs  and  structures  are  random.  This  is 
especially  true  for  protective  structures  in  a  blast  and  shock  environment. 
Dynamic  snock  inputs  are  often  considered  random  because  (i)  the  energy 
expended  at  the  shock  source  cannot  be  directly  measured,  and/or  (2)  the  phy¬ 
sical  properties  of  the  medium  connecting  the  snock  source  to  the  point  of 
interest  cannot  be  observed.  Structural  systems  are  considered  random  because 


i 


(1)  tne  parameters  in  tne  constitutive  laws  governing  behavior  of  the  struc¬ 
tural  materials  are  random,  and/or  (2)  the  structural  geometry  does  not 
precisely  match  the  intended  configuration.  (There  is  another  type  of  uncer¬ 
tainty  that  has  been  mistaken  for  randomness  and  that  is  modeling  accuracy. 
Moael  inadequacy  is  a  systematic  error  and  will  not  be  directly  addressed  in 
this  study.) 

The  purpose  of  a  reliability  analysis  is  to  compute  the  probability  of 
survival  of  a  structure.  Reliability  computation  is  difficult  when  the  struc¬ 
tural  response  is  nonlinear  and  has  large  deflections.  Therefore,  one  of  tne 
ultimate  goals  in  reliability  analysis  is  to  specify  a  technique  for  the  com¬ 
putation  of  reliability  of  structures  executing  nonlinear,  large  deflection 
response. 

The  purpose  of  this  investigation  is  to  establish  techniques  for  tne 
analysis  of  reliability  of  dynamically  loaded  structures  that  behave  in  a 
hignly  nonlinear  way.  A  progression  of  problem  complexity  will  cover  a  range 
of  techniques  from  single-degree-of -freedom  (SDF)  systems  subjected  to 
stationary  inputs  to  multiple  of-freedom  (MDF)  systems  subjected  to 

transient  loads.  There  are  four  requirements  to  these  approaches.  First,  the 
random  cnaracter  of  tne  structural  system  must  oe  identified.  Second,  tne 
random  character  of  tne  loads  must  be  described.  Third,  the  random  character 
of  tne  failure  criteria  must  be  established.  Finally,  the  fourth  requirement 
is  to  estaDlish  a  method  using  the  first  three  requirements  to  produce  a 
prediction  of  probability  of  failure  of  a  structural  system.  Necessar  ily, 
this  fourth  step  uses  an  intermediate  step  of  determining  the  probabi 1 istic 
nature  of  tne  response  and  then  applies  tne  random  failure  criteria  to 
establish  the  probability  of  failure  or  its  complement  reliability. 
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2.  LITERATURE  REVIEW 

There  are  several  areas  of  research  Interest  which  have  a  direct  bearing 
on  the  topic  under  consideration,  damage  analysis  of  randomly  excited  systems 
with  random  characteristics.  The  purpose  of  this  literature  review  is  to 
summarize  the  results  of  some  recent  studies  which  relate  to  the  present 
topic.  This  literature  review  deals  with  investigations  in  the  probabilistic 
theory  of  structural  dynamics,  and  investigations  that  characterize  the  appro¬ 
priate  loading  and  system  parameters  using  probabilistic  models.  Specific¬ 
ally,  first,  some  texts  and  papers  discussing  the  characterization  of  dynamic 
random  structural  response  will  be  listed.  Next,  specific  studies  useful  in 
the  analysis  of  failure  due  to  peak  response  are  given.  The  characterization 
of  damage  is  important  for  structures  which  can  fail  due  to  damage  accumula¬ 
tion;  therefore,  papers  in  this  area  are  discussed  next.  Then  the  studies 
which  take  advantage  of  the  Markov  property  of  a  structural  response  and  its 
measures  will  be  considered.  Finally,  some  papers  which  discuss  the  random 
character  of  environments  and  structural  systems  are  given. 

There  are  several  texts  which  deal  with  many  aspects  of  the  probabilistic 
theory  of  structural  dynamics.  Among  these  are  the  books  by  Lin  (Ref.  1), 
Cranoall  and  Mark  (Ref.  2),  Newland  (Ref.  3),  Crandall  (Refs.  4  and  5).  ana 
Clougn  ano  Penzien  (Ref.  6).  These  books  treat,  in  detail,  the  problem  of 
computation  of  response  moments  for  structural  systems.  -In  particular,  these 
booKs  present  techniques  for  finding  the  moments  of  the  structural  response  of 
linear  SOF,  MDF,  and  continuous  systems,  excited  by  stationary  and  nonstation¬ 
ary  random  inputs.  The  systems  considered  all  have  deterministic,  constant 
parameters.  Some  relatively  simple,  nonlinear  problems  are  also  considered. 


For  example,  some  approximate  means  for  computing  the  response  moments  of 
nonlinear,  elastic,  SDF  structures  are  given  in  Reference  1.  These  texts  also 
consider  such  problems  as  the  first  passage  problem  for  linear  SDF  systems, 
the  fatigue  problem,  and  the  Markov  character  of  the  response  of  systems 
excited  oy  wnite  noise. 

Many  general  papers  broadly  characterizing  one  or  more  aspects  of  struc¬ 
tural  response  to  a  random  input  have  been  written.  For  example,  the  paper  by 
Rice  (Ref.  7)  treats  the  response  moments  of  linear  and  simple  nonlinear 
systems.  It  also  characterizes  tne  frequency  of  passage  of  a  response  random 
process  beyond  a  barrier  of  fixed  height.  The  paper  by  Ang  (Ref.  8)  gives  a 
methoa  for  computing  the  response  moments  of  a  linear  system,  and  then  shows 
now  to  use  these  for  first  passage  ano  other  computations.  Tne  papers  oy 
Wirscning  and  Yao  (Refs.  9  and  10)  and  Bogdanoff,  Goldberg,  and  Bernard 
(Ref.  11)  evaluate  the  response  of  linear  structures  to  earthquake-type  dis¬ 
turbances  and  assess  structural  safety  in  light  of  the  response.  Uhlenbeck 
and  Urnstein  (Ref.  12)  and  Wang  and  Uhlenbeck  (Ref.  13)  consider  the  Markov 
character  of  the  response  of  linear  systems  to  white  noise  input  and  use  this 
to  derive  the  Foxker  Planck  equations  for  these  systems.  Tney  then  solve  the 
t-okker  Planck  equations  to  ootain  the  transition  probabi  1  ities  for  the  struc¬ 
tural  response.  Caugney  (Ref.  14)  derives  and  solves  the  Fokker  Planck  equa¬ 
tion  for  nonlinear  elastic  systems  and  obtains  the  response  moments  and 
transition  proDaoi 1 ities  for  the  structural  response.  Goldberg,  Bogdanoff, 
and  Sharpe  (Ref.  15)  and  Toland  and  Yang  (Ref.  16)  analyze  the  response  of 
simple,  nonlinear,  elastic  structures  and  evaluate  response  moments  and  first 
passage  probabilities.  Vanmarcke,  Yanev,  and  Oe  Estrada  (Ref.  17),  Paez  and 
Yao  (Ref.  la),  and  Iyengar  and  Iyengar  (Ref.  19)  characterize  tne  response  of 


SDF  hysteretic  structures.  They  compute  (1)  moments  and  first  passage 
probabilities  for  the  displacement  response,  and  (2)  the  moment  of  accumulated 
plastic  strain  and  permanent  set  in  the  structures,  hen  (Kef.  20)  specifies  a 
means  for  analyzing  MUF  hysteretic  structures.  He  places  each  hysteretic 
spring  witn  an  equivalent,  higher  order,  linear  spring,  kesponse  moments  are 
then  analyzed. 

Investigators  interested  in  identifying  the  probability  of  failure  of  a 
structure  often  consider  the  first  passage  problem,  or  equivalent  peak 
response  problem,  many  of  these  investigations  consider  SDF  systems.  For 
example,  the  papers  by  Yang  (Refs.  21  and  22),  Yang  and  Shinozuka  (Ref.  23), 
koberts  (Refs.  24  and  25),  Lin  (Ref.  26),  and  Corotls,  Vanmarcke,  and  Cornell 
(Ref.  27)  use  various  techniques  to  find  the  first  passage  probability  for 
linear,  SDF  oscillators.  Paez  (Ref.  28)  specifies  a  technique  for  computing 
tne  peak  response  prodaoility  distribution  for  an  SDF  system.  In  a  general 
paper,  Ang  (kef.  8)  defines  an  approximate  means  for  computing  the  first 
passage  probability  for  an  mUF  system.  Koopmans,  Qualls,  and  Yao  (Ref.  29) 
and  kojwithya  (kef.  30)  specify  bounas  on  the  peak  response  probability 
distribution  for  linear,  MQF  systems. 

In  assessing  the  chance  of  failure  of  a  structure,  the  potential  for 
accumulation  of  damage  must  oe  considered.  Several  authors  have  suggested 
codes  and  formulas  for  classifying  the  damage  in  a  structure.  For  example, 
qualitative  measures  of  structural  damage  are  discussed  in  papers  by  Whitman 
et  al.  (Refs.  31  and  32),  Housner  and  Jennings  (Ref.  33),  Hart  (Ref.  34),  and 
Hsu  (Ref.  35).  Quantitative  measures  of  the  damage  accumulated  in  simple 
structural  members  are  also  available.  For  example,  the  papers  by  Yao  and 
nunse  (kef.  36),  Tang  and  Yao  (Ref.  37),  Oliveira  (Ref.  38),  and  Kasiraj  and 
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Yao  (Kef.  39)  suggest  quantitative  measures  of  accumulated  structural  damage. 
Some  of  these  measures  of  damage  are  based  on  Palmgren-Mlner-type  damage 
accumulation  laws.  Some  literature  surveys  on  the  subjects  of  damage  assess¬ 
ment  and  measures  of  damage  have  been  written  by  Yao  (Refs.  40  and  41). 

in  investigations  into  the  probabilistic  structural  response  of  mech¬ 
anical  systems,  researchers  have  sometimes  taken  advantage  of  the  Markov 
character  of  the  response  to  white  noise  and  filtered  white  noise  Inputs, 
notably,  in  studies  of  the  first  passage  probability  for  linear  system 
response,  Yang  and  Shinozuka  (Refs.  42  and  43),  Rosenblueth  and  Bustamante 
(Kef.  44),  Gray  (Ref.  45),  Crandall,  Chandiramani,  and  Cooke  (Ref.  46),  and 
Paez  and  Yao  (Ref.  47)  nave  used  the  Markov  character  of  the  response.  All 
these  studies  concern  linear,  SOF  structures.  In  addition,  dogdanoff 
(Kefs.  48  and  49),  bogdanoff  and  Krieger  (Ref.  50),  and  Paez,  Tang,  and  Yao 
(Kef.  51)  have  showi  how  a  Markov  chain  approach  can  be  used  In  the  failure 
analysis  of  structural  systems.  When  the  damage  transition  probability  matrix 
is  provided,  and  when  the  probability  distribution  of  the  failure-causing 
damage  level  is  known,  their  approach  can  be  used  to  find  the  probability  of 
failure  of  a  structural  system. 

The  literature  review  in  this  section  summarizes  investigations  into 
various  aspects  of  the  probabilistic  structural  response  problem.  Many  papers 
dealing  witn  tne  probabilistic  theory  of  linear  structural  response  are  avail¬ 
able,  and  some  of  these  are  mentioned  in  review.  Fewer  papers  analyzing  tne 
peax  response  of  linear  structures  have  been  written,  and  even  fewer  deal  with 
nonlinear  structures.  Also,  only  a  few  papers  dealing  with  the  accumulation 
of  damage  in  randomly  excited  structures  are  available,  though  much  effort  is 
being  devoted  to  the  damage  accumulation  problem  at  this  time. 


14 


Damage  accumulation  models  may  also  be  affected  by  the  development  of 
more  sophisticated  constitutive  models.  There  is  extensive  work  being  done  in 
this  area  of  constitutive  modeling  of  concrete  soils.  Generally  the  failure 
criteria  will  be  more  adequately  specified  as  material  models  are  developed 
that  are  more  accurate  In  failure  regimes.  None  of  the  probabilistic  papers 
reviewed  considers  the  large  displacement  response  of  structural  systems.  A 
few  of  the  papers  in  the  literature  take  advantage  of  the  Markov  cnaracter  of 
structural  response;  most  of  these  consider  the  linear  response  behavior  of 
linear  systems.  This  property  of  structural  response  is  important  in  the 
development  of  this  investigation  and  can  be  extended  to  systems  that  respond 
in  a  nonlinear  way. 

The  random  cnaracter  of  tne  environment  for  pwotective  construction  has 
been  addressed  In  some  detail,  but  the  Information  is  not  available  in  the 
open  literature.  Finally,  there  are  some  papers  that  report  the  specific 
random  cnaracter  of  some  of  the  environment  and  structural  system  parameters. 
For  example.  Reference  52  summarizes  recent  studies  in  the  strength  of 


reinforced  concrete. 


II.  STOCHASTIC  ANALYSIS  OF  SIMPLE  MECHANICAL  SYSTEMS 


1.  INTRODUCTION 

In  this  section,  some  SDF  systems  are  analyzed  using  conditional  prob¬ 
ability  approaches.  Linear  elastic,  nonlinear  elastic,  and  hysteretic  SDF 
systems  are  analyzed.  The  parameters  of  most  of  the  systems  considered  are 
deterministic  constants;  however,  one  of  the  systems  analyzed  is  assumed  to 
have  a  random  capacity  for  damage  accumulation.  The  input  used  to  excite  most 
of  the  systems  considered  here  is  a  band-limited  white  noise  random  process. 

In  one  instance,  nowever,  a  nonstationary  random  process  is  used  to  excite  the 
system  under  consideration. 

Several  types  of  probaoilistic  results  are  obtained  for  the  SDF  systems 
analyzed  in  this  section.  First  passage  probabil itles  for  some  measures  of 
SuF  system  response  are  computed.  Finally,  the  failure  probability  is  com¬ 
puted  for  a  system  assumed  to  degrade  following  a  particular  damage  law. 

2.  FIRST  PASSAGE  OF  LINEAR  SDF  SYSTEMS 

a.  Theoretical  analysis— In  this  section  a  class  of  first  passage  prob¬ 
lems  for  SDF  systems  is  solved.  In  order  to  solve  the  first  passage  problem, 
tne  probability  that  some  measure  of  system  response  will  pass  outside  a  pre- 
astaolisned  barrier  at,  or  oefore,  a  specific  time  will  be  computed.  When  the 
response  level  corresponding  to  failure  in  an  SDF  system  can  be  defined  in 
terms  of  one  response  measure  and  when  it  is  deterministic  and  known,  tne 
probability  of  failure  of  an  SOF  system  can  be  computed  using  a  first  passage 
analysis. 

The  only  requirements  placed  upon  the  input  in  the  following  analysis 
are:  (1)  that  it  be  accurately  characterized  in  discrete  time,  and  (2)  that 


the  values  of  the  Input  random  process  at  consecutive  time  points  be  indepen¬ 
dent.  A  temporarily  stationary  Input  random  process  implies  that  the  input  is 
a  band-limited  white  noise.  In  the  analysis  of  SDF  systems,  the  assumption  of 
a  white  noise  Input  is  not  severe.  This  is  true  since,  in  general,  only  that 
value  of  the  spectral  density  of  the  input  near  the  natural  frequency  of  tne 
system  under  consideration  is  Important  in  Its  influence  on  the  structural 
response. 

first  passage  analyses  can  be  executed  using  a  Markov  chain  framework. 

The  only  quantities  that  must  be  generated  are  the  transition  probabilities. 
One  means  for  generating  transition  probaoilitles  will  be  presented  in  this 
section. 

The  equation  of  /notion  governing  tne  response  of  a  linear,  base-excitea 
SDF  system  is 


y  +  2>n  y  +  uij  y  *  -x 


(1) 


where  y  is  the  relative  displacement  response  of  tne  SOF  system,  x  is  tne  base 


acceleration,  <i»n  is  the  system  natural  frequency,  C  is  the  damping  factor,  and 


dots  denote  differentiation  witn  respect  to  time. 

In  order  to  obtain  the  transition  probabilities  required  for  the  proba¬ 
bilistic  analysis,  this  equation  is  first  discretized  in  time  using  a  central 
finite  difference  approximation.  The  resulting  approximation  to  the  governing 
equation  is 

:  a  *  y±i  + 
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Here  tne  equation  of  motion  is  written  for  time  t^  «  jAt,  j  *  0,  ...,N;  At  is 


tne  discrete  time  Interval  for  the  problem.  Equation  2  can  be  solved  for  y.+^ 
and  the  result  is 


(3) 
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When  the  system  parameters  are  constants  and  input ‘excitation  is  a 
stochastic  process,  the  SOF  system  response  is  also  a  stochastic  process. 

In  this  case  the  system  response  is  still  governed  by  Equation  3.  When 
the  actual  value  of  the  input  at  time  t.  is  specified  and  the  response 

J 

values  at  tj_^  and  tj  are  known,  the  the  response  value  at  t^+1  can  be 
approximately  computed.  In  the  following  it  has  been  assumed  that  the  system 
parameters  are  constants  and  the  input  excitation  is  a  stochastic  process. 

It  is  assumed  that  the  random  variables  composing  the  stochastic  process  are 
independent  at  consecutive  times,  and  that  the  probability  distribution  of  the 
input  is  known  at  each  time  tj ,  j  ■  0, . . .  ,N.  To  be  consistent  with  standard 
random  process  notation,  the  synbols  denoting  response  and  input  in  Equation  3 
are  capitalized  when  they  represent  stochastic  processes.  The  input  stochas¬ 
tic  process  is  denoted  Xj,  j  ■  0,...,  N;  the  response  stochastic  process  is 
Yj,  J  *  0,...,N. 

The  transitional  probabilities  governing  the  SDF  structure  response  will  be 
obtained  using  the  definition  of  conditional  probability  and  Equation  3. 

The  probability  that  the  structural  response  falls  in  the  interval  (y  ,  y.  ) 

a  0 

at  time  t^,  given  that  the  response  at  time  t^_1  is  y^  and  the  response  at 

time  t4  is  y,,  ,  is  given  by 
J  *i 

”(y* <  Vi  ±  yb|Vi  ‘  V  ’j  ’  yki) 

*  p(ya  <  (1  +  cwn  At)"1  [At2  Xj  +  (At2  U>2  -  2)  Yj  ♦  (1  -  ;un  At)  Y^] 
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The  expression  on  the  left  Is  simply  the  conditional  probability  described 
above.  The  expression  on  the  right  simply  uses  the  definition  for  Yj+j, 
established  In  Equation  3,  In  place  of  yj+j* 

Since  specific  values  for  the  response  at  times*  tj_i  and  tj  are  given 
in  Equation  4,  these  values  can  be  used  Inside  the  inequality  on  the  right- 
hand  side  of  Equation  4.  Moreover,  the  Inequality  can  be  solved  for  the 
Input  random  variable,  Xj.  The  resulting  expression  is 


P(ya  <  Vl  i  yblYj-l  '  yy  VJ  '  \) 
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When  specific  values  for  ya,  y^,,  y^g  and  y^  are  provided,  along  with  values 
for  the  system  parameters,  then  a  numerical  value  can  be  evaluated  for  the 
probability  In  Equation  5.  This  Is  true  since  the  probability  distribution 
of  Xj  has  been  assumed  known. 

To  establish  the  framework  for  the  computation  of  an  entire  collection  of 
probabilities  characterizing  the  structural  response,  the  present  problem  is 
now  discretized  in  response  space.  Let  -c  and  c^  represent  the  lower  and 
upper  first  passage  barriers.  The  first  passage  probability  Is  the  probabil¬ 
ity  that  Y.,  j  ■  0, . . .N  assumes  a  value  outside  the  Interval  (-c  ,  c.)  at  or 

J  «  D 

before  the  time  t^.  The  displacements  are  discretized  into  M  equally  spaced 
values  In  the  range  (-c  .  c.  ).  The  displacements  are 
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where 


yk  »  (k  -  1/2)  Ay  -  ca  , * 


Ay  »  (ca  +  cb)/M 


1.....M 
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(6a) 


Now  the  probability  that  the  response  at  tj+j  assumes  a  value  in  an 
interval  of  width  Ay  surrounding  the  displacement  yy  can  be  computed  by 
using  yj(  -  Ay/2  in  place  of  ya  and  y^  +  Ay/2  in  place  of  yb  in  Equa¬ 
tion  5.  When  the  response  at  tj+j  is  assuned  to  be  in  (y^  -  Ay/2,  yk 
+  Ay/2],  then  the  response  equals  y^.  This  assumption  makes  the  random 
variables  in  the  stochastic  process,  Yj,  j  *  0,...,N,  discrete  valued.  The 
assumption  results  in  an  accurate  representation  of  the  response  when  ay  is 
small  enough.  Mathematically,  the  assumption  described  above  is  written 

p(vi  ■  Vm  *  V  YJ  ■  ,kJ 

p(yk?  -  Ay n  <  Vl  <  •  %•  YJ  *  \) 

(7) 

The  numerical  values  of  all  the  transition  probabilities  represented  by 
Equation  7  can  be  computed.  In  all,  there  are  M3  values.  These  probabili¬ 
ties  characterize  the  transitions  of  the  response  from  points  within  the 

% 

first  passage  bounds  at  times  tj_j  and  tj  to  points  within  the  first  pas¬ 
sage  bounds  at  time  tj+j.  These  probabilities  do  not  represent  all  the 
possible  displacement  versus  time  paths  that  the  SDF  system  might  execute 
during  response  to  a  random  input.  Only  those  paths  remaining  within  the 
Interval  (-ca,  Cb)  Are  represented.  Those  transition  probabilities  which 
correspond  to  paths  which  originate  or  pass  outside  the  first  passage  bar¬ 
riers  will  not  be  used  in  the  present  analysis.  To  completely  characterize 
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Its  probabilistic  response,  starting  probabilities  for  the  SDF  system  under 
consideration  must  be  specified.  These  starting  probabilities  are  the  values 
of  the  joint  probability  mass  function  (pmf)  of  the  random  variables  Y. t  and 
Y0.  Y0  is  the  first  random  variable  in  the  response  stochastic  process;  Y_ l 
Is  a  random  variable  defined  at  the  ariflcial  time  t. l  *  -at.  Y_r  is  defined 
solely  to  facilitate  characterization  of  the  probability  distribution  of  the 
starting  velocity.  The  joint  pmf  of  Y_i  and  Y0  defines 


^0*^1  3  1 , . . . , M 


(8) 


Generally,  the  system  considered  will  have  a  zero  start  condition.  This 
means  that  the  system  starts  with  zero  displacement  and  zero  velocity.  This 
zero  start  condition  Is  characterized  by  the  joint  pmf 
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In  situations  where  none  of  the  discrete  displacements  actually  equals  zero, 
the  joint  pmf  is  set  to  one  at  the  argument  (y^,  yki)  where  y^  =  y^  is  the 
realization  of  the  random  variable  which  is  smallest  in  absolute  value. 

With  the  information  provided  In  Equations  5,  7,  and  3  or  9,  the  proba¬ 
bility  mass  function  of  the  response  stochastic  process  can  be  propagated 
through  time.  First  using  Equations  7  and  8  the  joint  pmf  of  Y0  and-Yl  can  be 
obtained  by  recognizing  that  multiplying  Equation  7  by  Equation  3  will  yield 
the  joint  pmf  of  Y_lt  Y0,  and  Yt.  Summing  over  all  k0  will  eliminate  depen¬ 
dence  Y.|  and  yield  the  Joint  pmf  of  Y0  and  Yr  Then  this  result  is  used  with 
Equation  7  and  the  same  process  to  get  the  joint  pmf  of  Yt  and  Y2,  etc.  The 
result  is,  at  time  tj_^  and  tj,  we  have  the  joint  pmf  of  and  Yj.  But 
recall  that  the  transition  probabilities,  developed  above,  consider  only  those 

oaths  remaining  in  the  spatial  interval  (-<:  .  c  '. 
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So,  ;n  eact,  the  luantity 


P* ( Y j_ ^  »  y^,  Yj  *  y^  ),  obtained  in  the  manner  described  above,  defines  the 
triple  joint  probability  that  (1)  the  response  at  time  tj_j  equals  yk  ,  and 
(2)  the  response  at  t^  equals  yfc  ,  and  (3)  the  SDF  system  response  has  not 
passed  outside  the  interval  ( -c  ,  c.  )  at  or  before  the  time  t..  (The  star  (*) 

a  D  J 

superscript  is  included  on  P*(Yj_^,  Yj  *  y^  )  to  show  that  only  a  portion  of  a 
complete  joint  pmf  is  defined  by  the  expression.  The  probabilistic 
description  of  paths  passing  outside  the  barrier  is  not  obtained.) 

Summing  this  expression  over  all  k0  yields 
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(10) 


This  is  the  joint  probability  that  the  response  at  tj  equals  and  the 
system  response  has  not  passed  outside  the  interval  (-ca,  C5)  at  or 
before  t j . 

When  Equation  10  is  sunmed  over  all  klt  the  "no  passage"  probability  is 
obtained.  If  T ^  is  defined  as  the  random  variable  denoting  the  time  when 
first  passage  of  the  SDF  system  relative  displacement  response  outside  the 
interval  (-ca,  c^)  occurs,  then 


P(Tj  > 
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This  is  the  chance  that  no  passage  of  the  structural  response  outside  the 
spatial  interval  (-ca,  C5)  occurs  in  the  time  interval  (0,  tj].  The 
first  passage  event  is  complementary  to  the  event  Tj  >  tj; 


P<T1  <  lj)  a  1  *  p(Tl  >  tj)  ,  j  -  0 . N 
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This  is  the  first  passaae  probability  function,  the  desired  quantity. 


Note  that,  in  order  to  use  the  present  equations  for  the  failure  analysis 
of  a  structure,  the  structure  must  behave  linearly  (or  nearly  linearly)  up  to 
the  failure  point,  and  it  must  be  assumed  that  failure  occurs  due  to  peak 
response.  Moreover,  the  system  must  be  accurately  modeled  as  an  SDF  system. 

On  the  other  hand,  the  input  can  be  modeled  quite  accurately.  The  input 
can  be  modeled  as  a  nonstationary  random  process  since  the  transition  proba¬ 
bilities  defined  in  Equation  5  depend  on  the  distribution  of  x^,  and  this 
distribution  may  vary  with  j. 

The  analysis  described  above  is  limited  to  the  consideration  of  first’ 
passage  problems  where  the  barriers  are  constant  and  define  a  spatial  interval 
( -c  ,  c. ).  In  fact,  the  analysis  can  be  generalized  to  make  the  time  domain 
discretization  vary  with  j.  One  could  simply  redefine  the  barriers  at  every 
time  tj  and  derive  the  discrete  displacements  as  in  Equations  6  and  6a.  This 
generalization  would  permit  the  performance  of  a  probabilistic  analysis  of 
failure  in  systems  where  the  strength  properties  degrade. 

Since  only  those  paths  remaining  between  the  first  passage  barriers  are 
considered  in  the  present  analysis,  only  a  portion  of  the  probabilistic 
structural  response  is  characterized.  In  fact,  the  entire  spectrum  of  poten¬ 
tial  responses  could  be  accounted  for  by  inclusion  in  the  analysis  of  two 
extra  discrete  displacements;  one  of  these  would  be  above  the  upoer  ’imit, 
Cfj,  and  the  other  would  be  below  the  lower  limit,  -ca.  All  paths  passing 
outside  the  first  passaae  barriers  at  a  time  tj  would  be  lumoed  in  these 
two  discrete  displacements.  The  probability  that  a  path  originates  outside  a 
barrier  and  then  returns  to  a  displacement  value  inside  the  barriers  would  be 
set  to  zero.  Using  this  approach  guarantees  that  the  probability  mass  func¬ 
tion  for  displacement  at  a  time  tj  would  have  elements  which  sun  to  unity. 
The  orobability  that  the  system  response  occupies  the  displacement  value 
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above  the  upper  limit,  Cb.  or  the  displacement  value  below  the  lower 
limit,  -ca,  at  time  tj,  would  be  the  first  passage  probability. 

Finally,  note  that  the  probability  distribution  of  the  stochastic  process 
input  has  been  left  indefinite  in  Equation  5.  In  many  situations  the  input 


will  be  defined  as  a  normally  distributed  stochastic  process,  and  in  such 
cases  Equation  5  will  be  specialized  to 
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*(•)  is  the  emulative  distribution  function  of  a  standard  normal  random 
variable;  uj  and  aj  denote  the  mean  and  standard  deviation  of  the  input 
random  process  at  time  tj. 


Numerical  examples— A  computer  program  which  executes  the  analyses 
developed  in  the  previous  section  was  written.  A  listing  of  the  computer 
program  is  included  with  this  report  In  Appendix  B.  Some  nunerlcal  examples 
demonstrating  the  use  of  the  computer  program  are  summarized  in  this  section. 

The  computer  program  computes  the  first  passage  probability  function  for 
a  base-excited,  linear  SDF  system.  The  input  excitation  is  a  stochastic  pro¬ 
cess  composed  of  a  sequence  of  independent  normal  random  variables.  Each 
Input  random  variable  has  mean  zero  and  the  standard  deviation  can  decay 
exponentially  with  time.  Zero  start  initial  conditions  for  the  SDF  system 


are  assumed. 


One  important  feature  of  the  program  is  the  amount  of  storage  required 
for  the  discrete  functions  used  in  the  analysis.  When  an  M-division  spatial 
discretization  is  used  in  the  problem  solution,  M3  storage  locations  are 
required  to  hold  the  transition  probabilities  of  Equation  7.  In  addition, 

M3  storage  locations  are  required  to  hold  the  products  between  the  joint 
probabilities  and  the  transition  probabilities  before  the  sunmation  Is  exe¬ 
cuted.  (Improvements  in  the  analytical  approach  which  diminish  the  required 
amount  of  storage  are  described  later.)  Other  functions  increase  the  storage 
requirements  as  functions  of  M  and  M2;  therefore,  more  than  2M3  storage  loca¬ 
tions  are  needed  for  this  program.  Clearly,  this  limits  the  range  of  discret¬ 
izations  that  can  be  used  in  problem  solutions. 

The  program  listed  in  the  Appendix  was  actually  written  to  solve  time  and 

space  normalized  versions  of  the  problem  discussed  In  the  previous  section. 

In  particular,  the  program  performs  a  natural  frequency  independent  computa- 

* 

tion,  where  the  probability  of  first  passage  at  or  before  a  nunber  of  cycles 
of  SDF  system  response  is  determined.  In  space,  the  upper  and  lower  first 
passage  bounds  are  written  as  a  multiple  of  the  root  mean  sauare  (rms)  value 
of  the  stationary  response  of  an  SDF  system  when  the  input  is  stationary. 

The  first  problem  to  be  summarized  is  the  first  passage  of  a  base 
excited,  linear  SDF  system.  The  input  excitation  is  a  stationary  random 
process.  Specifically,  it  is  a  band-limited  white  noise.  The  standard  devia¬ 
tion  of  the  Input  is  listed  In  Table  1.  The  system  and  computation  param¬ 
eters  are  listed  in  Table  2.  The  first  passage  probability  was  determined 
at  each  computation  cycle  and  the  results  are  shown  in  Figure  1.  (Note 
that  curves  relating  to  two  examples  are  given  in  Figures  1  and  2.  The 


curves  marked  "stationary"  are  connected  witn  this  example;  curves  marked 
"nonstationary"  are  connected  to  a  later  example.)  The  first  passage  proba¬ 
bility  Is  graphed  as  a  function  of  time  In  Figure  1.  In  this  example,  the 
first  passage  probability,  at  a  particular  time,  Is  the  chance  that  the  dis¬ 
placement  response  of  the  structure  has  passed  outside  the  spatial  interval 
(-ca,  cb),  at  or  before  that  time.  As  expected,  the  first  passage  probability 
Increases  with  time.  Other  examples,  not  summarized  here,  show  that  when  the 
barrier  values  are  Increased  in  absolute  value,  the  first  passage  probability 
function  Increases  more  slowly  with  time.  When  the  barrier  values  are 
decreased,  the  first  passage  probability  Increases  more  rapidly  with  time. 

When  the  excitation  stochastic  process  is  stationary,  the  first  passaae 
probability  always  approaches  unity  asymptotically,  as  time  increases.  The 
reason  is  that  a  certain  portion  of  the  probability  for  response  paths  within 
the  barriers  Is  lost  to  paths  outside  the  barriers  at  each  step  in  the 
computation.  When  the  input  is  nonstationary  and  decaying,  this  is  not 
necessarily  so. 


The  second  problem  to  be  summarized  also  concerns  the  first  passage  of  a 
base-excited,  linear  SOF  system.  In  this  case,  however,  the  input  is  a  non¬ 
stationary  random  process.  The  input  is  a  sequence  of  Independent,  mean 
zero,  normally  distributed  random  variables  with  standard  deviations  that 
vary  as  a  function  of  time.  Specifically 
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This  function  is  graphed  in  Figure  7.  along  with  the  constant  value  of  stan¬ 
dard  deviation  from  the  previous  problem.  The  first  passage  bounds  were  cho¬ 
sen  equal  in  maanitude  to  those  used  in  the  previous  problem. 

The  first  passaqe  probability  is  graphed  as  a  function  of  time  in  Figure 
1.  Since  the  input  used  in  this  problem  is  less  severe  than  that  used  In  the 
previous  problem,  the  first  passage  probability  is  lower  at  all  times. 
The  first  passage  probability  computed  here  does  not  approach  the  value  unity 
as  t  Increases.  This  Implies  that  first  passage  in  this  system  response  is 
not  certain.  The  reason  for  this  behavior  is  that  the  input  becomes  negligi¬ 
ble  as  time  progresses. 


c.  Computational  enhancements — The  conditional  probabilities  that  dis¬ 
placements  yv2,  k 2  *  l,...,u  will  be  realized  at  time  tj+j,  given  the 
displacement  combinations  at  two  previous  time  steps,  y(<0  and  yi^,  k0,ki 
»  1, . . . ,M  at  times  tj.j  and  tj  can  be  assembled  In  a  computer  code  as  a 
three-dimensional  array.  When  displacements  are  discretized  into  M  spatial 
increments  then  the  array  becomes  an  M  by  M  by  M  array  requiring  M3  storage 


locations.  As  discussed  previously,  this  array  is  treated  as  a  tensor  and  is 
multiplied  by  the  state  probability  matrix  which  is  composed  of  the  probabil¬ 
ities  that  the  disol acements  actually  were  at  specific  discrete  locations  at 


the  two  previous  time  steps.  The  multiplication  yields  the  joint  probabili¬ 
ties  that  the  displacements  are  at  three  locations  at  times  t.  . ,  t.,  and  t...» 

j-i  j  j+i 

and  these  probabilities  are  also  stored  In  an  array  of  size  M3.  This  multi¬ 
plication  Is  conceptually  shown  In  Figure  3.  The  next  step  In  the  analytical 
process  Is  to  sum  the  vectors  In  the  j-1  direction  into  scalars.  This  pro¬ 
duces  an  M  by  M  array  that  Is  the  joint  probability  that  the  displacements  are 
at  discrete  locations  at  times  tj  and  Vl*  Note  that  these  form  the  state 
probability  matrix  for  the  next  Incremented  time  step  (where  j-*j-l  and  j+l+j). 
The  summing  process  Is  conceptually  shown  In  Figure  4. 

By  simply  noting  that  matrix  multiplication  is  a  series  of  vector  opera¬ 
tions  It  can  be  shown  that  only  a  single  vector  of  the  conditional  probabil¬ 
ity  tensor  need  be  formed  at  one  time.  The  shaded  portions  of  the  conceptual 
figure  illustrate  that  a  single  vector  of  the  conditional  probability  tensor 
multiplied  by  the  state  probability  matrix  yields  a  single  vector  of  .the 
joint  probability  tensor.  This  vector  can  be  sunmed  into  a  scalar  element  of 
the  state  probability  matrix  for  the  next  time  step  (Figure  4).  When  this 
ooservatlon  is  used  in  the  computer  coding  algorithm,  the  storage  requirements 
for  the  computational  process  are  suostantially  reduced.  In  particular,  the 
storage  requirement  is  reduced  from  2M3  +  M2  locations  to  M2  +  2M  locations. 

wnen  the  input  excitation  is  stationary,  the  conditional  proDaoility  ten¬ 
sor  is  composed  of  elements  that  are  constant  throughout  the  computation. 
Therefore,  the  approach  described  above  requires  the  computation  of  these 
constants  each  calculatlonal  time  step.  This  means  that  storage  requirements 
are  reduced  at  an  increase  In  computational  time.  However,  the  computations 
performed  in  thfs  project  were  done  on  a  CRAY  computer  and  the  elements  are 
formed  as  a  vector  operation  on  the  CRAY  comDuter.  therefore,  the  'ncreaseo 


computational  time  Is  not  severe.  Furthermore,  many  of  the  Inputs  which  are 
of  Interest  are  nonstationary  and  the  conditional  probability  matrix  must  be 
recomputed  each  time  step  regardless  of  the  scheme.  Therefore,  the  above 
change  in  the  calculational  process  represents  an  important  step  in  reducing 
storage  requirements.  In  fact,  this  general  approach  may  allow  the  process 
to  be  expanded  to  multiple  degrees  of  freedom  without  surpassing  computer 
storage  requirements. 

Other  computational  enhancements  are  available.  For  example.  It  would  be 
much  more  efficient  to  compute  first  passage  probabilities  If  single-step 
memory  were  required  rather  than  the  two-step  memory  shown.  Recall,  the  two- 
step  memory  came  about  as  a  result  of  the  central  difference  scheme  used  to 
approximate  the  differential  equation.  One  method  to  diminish  the  memory 
would  be  to  use  an  approximation  that  would  require  only  single-step  memory, 
but  the  accuracy  would  suffer.  Another  approach  Is  to  Identify  special  fea¬ 
tures  of  the  response  such  that  tne  memory  may  be  reduced.  An  example  is  that, 
for  lightly  damped  SOF  systems,  the  response  to  wnite  noise  Input  is  narrow 
oand.  That  is.  given  that  the  response  is  at  a  peak,  at  a  given  time,  the  time 
at  wnich  it  will  be  a  minimum  will  be  approximately  one-half  the  natural  period 
of  the  system  later.  This  half-cycle  response  information  can  be  used  to 
reduce  the  computational  procedure  to  single-step  memory. 

To  demonstrate  this  approach  the  scheme  developed  above  has  been  used  to 
compute  the  probability  distributions  of  response  amplitudes  after  a  half 
cycle  of  response,  given  that  the  start  conditions  are  zero  velocity  and  any 
of  the  discrete  displacements  defined  In  Equation  6.  These  probabilities  are 
known  as  half-cycle  transition  probabilities,  and  since  the  analysis  considers 
only  those  paths  remalnlnq  In  the  spatial  Interval  (-ca,  c^),  these  probabili¬ 
ties  have  special  meaning.  Specifically,  these  probabilities  define  the 
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chance  that  the  SOF  system  response  Mill  be  at  an  amplitude  value  a.  at  the 

*  1 

end  of  a  half-cycle  and  that  the  response  Mill  not  pass  outside  the  Interval 

(-c  ,  c^)  during  the  half  cycle,  given  that  the  response  started  Mlth  velocity 

zero  and  an  amplitude  a.  at  the  beginning  of  the  half-cycle. 

X0 

The  sequence  of  mathematical  operations  used  to  obtain  the  half- cycle 
transition  probabilities  is  as  follOMS.  Set  the  starting  probabil ities. 
Equations  8  and  9,  to 
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Next,  divide  a  response  half-cycle  (duration  into  an  integral  number 

of  parts,  n‘.  Apply  the  process  n‘  times,  starting  Mlth  pY  v  (y.  ,  y.  ) 
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above  to  get  PY  Y  (yv  yk  ).  Finally,  apply  Equation  10  to  obtain  the 
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distribution  of  response  values  after  a  half-cycle,  considering  only  those 
paths  Mhlch  remained  betMeen  the  barriers. 
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Let  Aj,  j  «  0, ...,N'  be  the  SDF  system  amplitude  response  stochastic 
process,  and  let  Aj  be  the  random  variable  in  the  process  denoting  the 
response  amplitude  at  half-cycle  j.  Then  the  half-cycle  transition  probabil¬ 
ities  derived  using  the  scheme  described  above  are 
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where  the  right-hand  side  1$  obtained  In  Equation  16.  Since  the  computation 

scheme  has  been  used  to  compute  the  transition  probabilities  given  above.  It 

Is  convenient  to  choose  the  values  of  a,.  and  ab  as  those  defined  In  Equation 

6.  The  star  {*)  superscript  Is  Included  on  the  expressions  written  above  to 

denote  the  fact  that  only  the  response  paths  remaining  In  the  spatial  Interval 

(-c.,  O  are  Included  in  the  computation.  Paths  which  connect  au  to  a,,  in 
a  o  *  i 

Equation  17,  but  pass  outside  (-ca,  c^)  during  the  half-cycle  are  neglected. 

The  transition  probabilities  presented  In  Equation  17  provide  the  s&me 
sort  of  information  as  provided  by  the  two-step  memory  transition  probabili¬ 
ties  In  Equation  7.  The  obvious  difference  Is  that  here  only  one  memory 
step  is  included.  The  need  for  one  item  of  Information  in  the  memory  has 
been  eliminated,  since  only  the  response  at  half-cycle  points  is  considered. 

At  both  the  beginning  and  end  of  every  half-cycle  the  velocity  is  zero; 
therefore,  the  velocity  values  need  not  be  considered.  When  a.  two-step 
memory  is  used,  both  (average)  velocity  and  displacement  at  times  tj_i  and 
tj  influence  the  response  at  tj+i*  The  disadvantage  in  considering  only 
the  response  at  half-cycle  points  is  that  a  half-cycle  duration  must  be 
chosen  for  use  in  the  computation  to  obtain  the  probabilities  in  Eauation 
17.  This  is  chosen  as  n/u^  for  the  SOF  system.  While  this  is  the  average 
half-cycle  duration  for  a  system  responding  to  a  random  input,  the  actual 
half-cycle  duration  is  random.  Some  error  is  committed  in  using  one  value 
for  the  half-cycle  duration,  but  when  damping  is  light  (c  <  0.20)  the  error 
is  small  and  is  more  than  justified  In  reducing  the  memory  requirements. 

Figure  5a  shows  the  distribution  of  response  amplitudes  after  a  half¬ 
cycle  of  response  when  the  upper  and  lower  limits  are  chosen  as  c^  »  4ov 
and  -ca  «  -Aoy,  and  the  starting  amplitude  is  chosen  as  3.8  ay.  Also 
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shown  is  a  displacement  path  the  response  miqht  execute  during  the  half 
cycle.  The  lines  on  the  right  represent  the  probabilities  that  the  various 
amplitudes  will  be  realized  and  the  response  path  will  remain  in  (-4oy, 
40y)»  given  that  the  starting  amplitude  is  3.8  o^.  These  probabilities  form 
a  vector  of  transition  probabilities  given  the  starting  amplitude  of  3.8  ay. 

Figures  5b  and  5c  show  other  distributions  of  amplitude  response  proba¬ 
bilities  given  starting  amplitudes  of  2.2  ay  and  0.2  oy.  These  graphs 
present  information  similar  to  that  described  above. 

Half-cycle  transition  probability  vectors  corresponding  to  each  of  the 
beginning  amplitudes  can  be  generated.  These  vectors  can  be  combined, 
sequentially,  based  on  tne  starting  amplitude,  and  the  result  is  a  transition 
probability  matrix.  Only  those  paths  remaining  in  the  spatial  interval  (-c  , 

u 

cb)  are  included.  The  transition  matrix  is  defined  as  follows. 
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The  values  of  the  pmf  for  the  various  amplitude  states  at  half-cycle  zero  form 
a  vector  {Pg},  where  (Pg}  is  defined  as 
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The  star  (*)  superscript  indicates  that  only  the  probability  of  those  ampli¬ 
tudes  within  the  interval  (-ca,  cfa)  are  included  in  {Pg};  therefore,  the 

★ 

sun  of  the  elements  in  {Pg}  may  not  be  one.  Then  the  distribution  of  ampli¬ 
tudes  at  half  cycle  j  is 
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where  only  those  paths  remaining  in  the  spatial  interval  (-ca,  c^  are 
considered  in  the  computation. 

The  probability  that  no  passage  of  the  SDF  system  response  outside  the 
interval  (-ca,  Cb)  occurs  between  half-cycle  zero  and  half-cycle  j  is 


PI 
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Tl  is  the  first  passage  time.  p.  (a.),  k  *  are  the  elements  in  the 

Aj  K 

vector  {p j } .  Equation  74  yields  the  no-passage  probability  for  the  time 
interval  (0,  *j/wn)*  In  terms  of  this  quantity,  the  first  passage  probability 
can  be  written 

P(T1  <  irj/«n)  »  1  -  P(T1  >  *j/»n)  ,  j  »  1 . N1  (22) 
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The  first  passage  probability  can  be  computed  as  a  function  of  time  when 
the  procedures  outlined  above  are  followed.  The  computations  performed  using 
these  half-cycle  transition  probabilities  are  subject  to  the  assumptions  made 
In  the  previous  section. 

When  the  Input  excitation  stochastic  process  is  stationary  in  time,  the 
present  technique  for  analysis  of  first  passage  probabilities  of  a  linear, 
base-excited  SDF  system  will  yield  a  computational  time  savings  proportional 
to  MN/N'.  This  factor  will  be  considerable  when  the  spatial  discretization 
is  fine  (M  is  large)  and  the  time  discretization  is  fine  (N/N*  is  large  and 
both  analyses  cover  the  same  interval  of  time).  When  the  input  is  nonsta¬ 
tionary,  the  present  technique  does  not  result  in  computational  time  savings. 

Expansion  of  the  above  analytic  approach  to  highly  nonlinear  or  MDF 
problems  is  not  direct.  Tne  oasic  assumption  of  narrow  band  response  would  not 
in  general  be  true  in  these  cases.  However,  the  concept  of  identifying  special 
features  of  the  response  or  the  computational  scneme  holds.  For  example, 
because  of  stability,  the  displacement  of  a  node  in  a  finite  element  scheme  is 
bounded  by  the  time  step.  There  are  certain  displacements  that  have  absolutely 
zero  probability  of  being  reacned  In  a  single  time  step.  Therefore,  it  may  be 
possible  to  extend  the  above  concept  into  solving  the  much  more  difficult 
problems  of  nonlinear  MOF  problems. 

3.  FIRST  PASSAGE  OF  NONLINEAR  ELASTIC  SO F  SYSTEMS 

a.  Theoretical  analysls--Once  the  first  passage  problem  for  a  linear 
elastic  SOF  system  has  been  solved,  it  Is  a  relatively  simple  matter  to 
extend  the  analysis  to  nonlinear  elastic  SOF  systems.  The  analysis  presented 
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above  Is  extended  in  this  section  to  include  the  nonlinear  elastic  case.  The 
input  used  to  excite  the  SDF  system  is  the  same  as  that  used  previously;  that 
ist  the  input  is  a  sequence  of  independent,  mean  zero  random  variables. 

The  equation  of  motion  governing  the  response  of  a  nonlinear  elastic, 
base-excited  SOF  system  can  be  written 


y  +  2ca»n  y  + 
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where  y  is  the  relative  displacement  response  of  the  SOF  system,  x  is  the 
base  acceleration,  is  the  "small  displacement  response"  natural  fre¬ 
quency  of  the  system,  c  is  the  damping  factor  and  dots  denote  differentiation 
with  respect  to  time.  This  equation  governs  the  response  of  a  system  whose 
restoring  force  is  an  odd  polynomial  function,  of  degree  2n  -  1. 

To  obtain  from  this  equation  the  transition  probabilities  required  for  a 
probabilistic  response  analysis,  first  discretize  the  equation  in  time.  When 
the  central  difference  approximations  are  used  to  replace  the  derivatives  in 
the  above  expression,  then 
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This  equation  governs  system  motion  at  times  tj  =  jAt,  j  *  0,...,N.  Re¬ 
arrangement  of  this  equation  yields  an  expression  for  yj+i. 


yj+l  *  t1  +  ?“n 


-1 


l“‘  5j  MK  “n  &  Vj 


+  (1  -  ;u»n  At) 
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Wien  a  nonlinear  elastic  system  is  excited  by  an  input  stochastic  pro¬ 
cess,  the  response  at  tj+^  is  approximately  governed  by  the  above  equation  if 
the  system  parameters  are  constants.  To  represent  the  approximate  stochastic 
relationship  between  input  and  response,  the  lower  case  yjS  must  be  replaced 

M  *  •• 

with  upper  case  Y.s  and  the  lower  case  x.  with  Xi.  Then  using  the  definition 

J  J  J 

of  conditional  probability,  an  expression  for  the  SDF  system  response 
transition  probabilities  can  be  obtained.  Tnat  is, 

p(^<  yii'Jvi-  vvy0 


kl  i 
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This  expression  is  the  nonlinear  system  equivalent  to  the  transition  proba¬ 
bility  expression  derived  in  Section  II. 2. a  for  linear  systems;  and  it  is 
derived  in  the  same  way  and  can  be  used  in  exactly  the  same  way.  It  is  noted 
that  not  only  stability  of  the  finite  difference  equation  must  be  properly 
addressed,  but  that  adequate  time  step  is  selected  to  propagate  probability 
outside  initial  conditions. 

When  the  system  response  displacements  are  discretized.  Equation  26  can 
be  used  to  obtain  the  transition  probabilities  for  the  discrete-time/discrete- 
space  SOF  system  response  stochastic  process.  The  assumptions  and  notations 
regarding  response  displacement  bounds  are  the  same  here  as  those  used  in 

Section  II. 2. a.  The  discrete  displacements  in  (-ca>  c^)  are  y1<?  k  *  l . M 

defined  in  Equation  5.  The  system  response  transition  probaDi Titles  are 


p(vi  ■  yk2|yj-i  ■  V  vj  ■  ,ki) 


•  pK  -  *  iy/2b-i  ■  %•  yj  ■  \Y  ;;kX':"i . m 


where  the  expression  on  the  right-hand  side  is  evaluated  using  Equation  26. 

Now  when  the  starting  joint  pmf  of  the  response  stochastic  process  is 
specified  as  in  Equation  8  or  9  and  the  probability  computation  is  marched  out 
In  time,  as  before,  first  passage  results  equivalent  to  those  obtained  In  ‘ 
Section  II. 2. a  are  obtained  for  the  present  problem.  The  analysis  required 
here  is  precisely  equivalent  to  that  described  in  Section  II. 2. a,  following 
Equation  9.  Some  nunerical  examples  are  presented  in  the  following  section  to 
demonstrate  the  results  of  this  analysis. 

In  order  to  use  the  present  approach  for  the  practical  failure  analysis 
of  a  real  system,  it  Is  necessary  that  (1)  the  system  be  accurately  modeled 
as  a  nonlinear  elastic,  SDF  system,  and  (2)  failure  occurs  due  to  peak 
response. 

The  computational  enhancements  discussed  in  Section  II. 2. c  can  be 
applied  in  the  nonlinear  first  passage  analysis.  The  first  technique  dis¬ 
cussed  in  that  section  can  be  applied  without  modification.  The  second  com¬ 
putation  scheme,  the  conversion  of  the  response  random  process  from  a  condi¬ 
tional  probability  random  process  with  two-step  memory  to  one  with  one-step 

memory,  can  be  applied  here,  but  only  in  a  restricted  way.  An  underlying 

assunptlon  in  the  establishment  of  that  computation  scheme  is  that  the 

response  will  be  narrow  band.  When  the  degree  of  nonlinearity  in  SDF  system 
response  Is  not  too  great,  and  when  system  damping  is  light,  then  the 

response  will  be  narrow  band.  When  the  response  nonlinearity  becomes  great, 
the  response  will  cease  to  be  narrow  band.  Ductility  ratios  of  greater  than 
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1.5  and  damping  coefficients  of  greater  than  20  percent  of  critical  would 
be  considered  limits  for  these  approaches. 

The  following  analyses  characterizing  the  response  of  nonlinear  SDF  sys¬ 
tems  will  make  the  assumption  of  narrow  band  response  and  use  the  computa- 
ional  procedures  previously  discussed. 

b.  Numerical  examples— The  first  computer  program  listed  in  the  appen¬ 
dix,  and  used  to  solve  the  linear  first  passage  problems,  can  also  be  used  to 
solve  nonlinear  elastic  first  passaqe  problems.  Some  numerical  examples 
demonstrating  the  capability  of  the  computer  program  to  solve  nonlinear  first 
passaqe  problems  are  presented  in  this  section. 

The  computer  proqram  computes  the  probability  of  first  passage  of  the 
response  of  a  nonlinear  elastic,  base-excited,  SDF  system,  outside  the  spa¬ 
tial  interval  (-ca,  cb) .  The  response  is  excited  by  a  stochastic  process 
input  composed  of  a  sequence  of  zero  mean,  normal  random  variables.  Zero 
start  initial  conditions  for  the  SOF  system  are  assumed. 

Three  first  passaqe  problems  are  presented  in  this  section.  The  system 
parameters  used  in  Equation  23  through  27  are  sunmarized  in  Table  2.  The 
input  is  a  band-limited  white  noise  whose  standard  deviation  is  qiven  in 
Table  2.  The  first  passaqe  probability  was  computed  for  each  of  the  three 
systems  (linear,  hardening,  and  softening),  and  the  results  are  shown  in 
Figure  6. 


TABLE  2.  SYSTEM,  INPUT  AND  COMPUTATION  PARAMETERS 
FOR  NONLINEAR  FIRST  PASSAGE  PROBLEMS 
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Linear  system:  m  »  1,  cj  * 

1 
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Nonlinear  hardening  system: 

m  *  2,  Ci  *  1,  C2  ■  1.5 
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Nonlinear  softening  system: 

m  «  2,  Ci  ■  1,  C2  *  1.5 

0.714 

0  0.5  1.0  1.5  2.0  2.5  3. 


Time  (seconds) 

Figure  6.  First  passage  probability  for  linear,  hardening  and 
softening  SDF  systems  wn  *  2*. 


The  hardening  system  tends  to  have  lower  response  amplitudes  than  the 
linear  system;  therefore,  the  first  passaqe  probability  curve  of  the  harden¬ 
ing  system  is  below  the  first  passage  probability  curve  of  the  linear  sys¬ 
tem.  The  softening  system  tends  to  exhibit  higher  response  amplitudes  than 
the  linear  system;  therefore,  the  first  passaqe  probability  curve  of  the 
softening  system  falls  above  the  first  passaqe  probability  curve  of  the 
linear  system. 

4.  PROBABILITY  DISTRIBUTION  OF  SYSTEM  RESPONSE 

a.  Theoretical  analysis — Up  to  this  point  only  the  first  passage  proba¬ 
bility  has  been  used  to  characterize  structural  response.  The  first  passage 
problem  has  been  oursued  since  (1)  peak  response  can  be  an  important  crite¬ 
rion  in  predicting  system  failure,  and  (2)  relatively  little  work  has  been 
done  toward  solving  first  passage  problems.  On  the  other  hand,  much  work  has 
been  done  toward  characterizing  the  probability  distribution  of  the  overall 
response  of  a  linear  system.  This  section  modifies  the  general  approach 
outlined  In  Sections  II. 2  and  II. 3  to  compute  the  approximate  probability 
dlstrioution  of  the  overall  response,  rather  than  just  the  first  passage 
probability. 

The  main  reason  it  is  difficult  to  comDute  the  probability  distribution 
of  the  overall  structural  response  is  that  the  procedures  outlined  in  the 
previous  sections  are  ideally  suited  to  tracking  a  portion  of  the  response 
probability  remaining  between  finite  barriers.  The  probability  connected 
with  paths  passing  outside  the  barriers  is  discarded  in  succeeding  computa¬ 
tions.  When  the  probability  distribution  of  the  overall  response  is  of 
Interest,  all  paths  must  be  accounted  for  and  none  of  the  probability  can  be 
discarded,  no  matter  now  great  the  amplitudes  of  the  paths  involved. 


Since  the  computational  scheme  used  with  the  present  approach  cannot 
account  for  all  paths,  especially  those  whose  amplitudes  are  very  large,  a 


compromise  Is  made  by  taking  two  actions:  (1)  by  setting  the  barriers  of  the 
spatial  Interval  (-c  ,  c.  )  at  very  large  amplitudes  so  that  the  probabilities 

a  D 

associated  with  paths  crossing  outside  the  barriers  are  low,  and  ( Z )  by  normal¬ 
izing  the  probabilities  to  a  value  of  1  at  every  step  of  the  computation; 

that  is,  by  changing  the  probabilities,  p*  v  (y,,  ,  yu  ),  referred  to  in 

j-1  j  11  o 

Equation  10.  On  the  jth_ computation  Interval,  the  sum 
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Is  the  no-passage  probability  through  computation  cycle  j. 

*  *  -1 
PY  v  *  yv  )  is  modeled  by  multiplying  it  by  (S  )  , 

Tj-1  j  0  1 
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If 

the  result  is 
kx  =  1 . M 


(29) 

The  path  probabilities  have  been  modified  so  that  the  chance  that  the 
reSDonse  remains  in  (-ca>  C5)  is  one.  When  the  interval  (-c3,  C5)  is 
wide  enouqh,  the  probability  distribution  of  the  response  will  be  accurate. 
The  star  (*)  has  been  dropped  on  the  left,  above,  since  now  the  probability 
for  all  response  paths  has  been  accounted  for,  thouoh  in  an  approximate  way. 

The  distribution  of  the  response  can  now  be  computed  at  each  time  point, 
and  is  obtained  by  summing  the  expression  in  Equation  29  over  all  k0. 


(30) 


46 


The  value  of  the  above  modification  in  computation  scheme  does  not  lie  in 
the  capability  it  establishes  to  approximate  the  distribution  of  the  response 
of  a  linear  SOF  system.  This  capability  already  exists  using  a  theoretical 
probability  analysis  of  linear  continuous  systems.  Rather,  its  value  lies  in 
the  fact  that  the  present  analysis  can  be  extended  to  nonlinear  elastic  and 
hysteretic  systems.  This  extension  will  be  outlined  in  following  sections. 
The  fact  that  the  method  presented  above  yields  results  corresponding  to  an 
exact  analysis  will  be  used  to  check  the  accuracy  of  the  present  approach. 

b.  Numerical  examp1e--The  analysis  outlined  in  Section  II. 4. a  is 
checked  for  accuracy  in  this  mmerical  example.  Three  numerical  examples 
were  solved.  In  each  example  the  response  probability  distribution  was  com¬ 
puted.  The  input  is  a  band-limited  white  noise  random  process  whose  standard 
deviation  is  given  in  Table  3.  After  tne  response  reacned  a  steaay  state,  the 
response  displacement  pmf  was  recorded.  The  parameters  of  the  tnree  problems 
solved  are  given  in  Table  3.  In  each  problem  tne  lower  and  upper  spatial 
limits  are  defined  in  terms  of  rms  response,  once  it  nas  reached  the  stationary 
state. 

Figure  7  shows  plots  of  the  response  cdf's  for  each  problem.  Also  shown 
is  the  theoretical  cdf  for  the  linear  SDF  structural  response.  The  statis¬ 
tics  of  the  theoretical  response  are  available,  for  example,  in  Reference  1. 
The  plots  show  that  there  is  reasonable  agreement  between  the  theoretical 
results  and  the  results  obtained  in  Section  II. 4. a. 

5.  DISTRIBUTION  OF  OISSIPATEO  ENERGY  IN  AN  ELASTO-PIASTIC  SYSTEM 

a.  Theoretical  analysis— The  analysis  of  structural  reliability  hinges 
on  our  ability  to  identify  the  modes  of  failure  which  might  lead  to  the 


TABLE  3.  SYSTEM,  INPUT  AND  COMPUTATION  PARAMETERS  FOR  RESPONSE 
PROBABILITY  DISTRIBUTION  PROBLEMS 


demise  of  a  structural  system.  Modes  of  structural  failure  can  be  qualita¬ 
tively  identified,  and,  to  some  extent,  quantitative  predictions  of  failure 
can  be  made.  8ut  to  date,  the  time  of  failure  cannot  be  exactly  predicted, 
even  when  the  response  history  is  computed  exactly.  The  reason  is  that  mate¬ 
rial  behavior  is  random.  Yet  in  a  reliability  analysis  some  criterion  must 
be  chosen  to  judge  whether  or  not  failure  occurs  in  a  response.  The  charac¬ 
terization  of  the  criterion  can  be  deterministic  or  probabilistic. 

Consider  the  el asto-pl ast ic  SDF  system.  Various  criteria  can  be  used  to 
judge  failure  in  this  system.  For  example,  accunulated  plastic  deformation 
can  be  used  to  determine  when  failure  will  occur,  or  a  permanent  set  can  be 
used  as  a  failure  criterion.  The  amount  of  enerqy  dissipated  by  the  system 
can  be  used  as  a  failure  criterion,  or  a  combination  of  these  factors  can  be 
used  as  a  failure  criterion.  In  the  present  analysis,  the  amount  of  energy 
dissipated  by  the  system  is  considered.  There  are  two  reasons  for  this:  (1) 
the  analysis  of  probability  distribution  and  first  passage  probabilities  for 
this  quantity  is  slightly  more  difficult  than  the  equivalent  analysis  for  the 
other  quantities;  and  (2)  this  might  serve  as  a  damage  criterion  for  an 
actual  system. 

The  equation  governing  motion  of  a  base  -  exc ited,  el asto-pl ast  ic ,  SDF 
system  is 

y  +  2qwn  y  +  u»^  R(y)  =  -x  (21) 

where  y  is  the  relative  displacement  of  the  system  mass,  x  is  the  base  accel¬ 
eration,  ay,  is  the  small  displacement  natural  frequency  of  the  system,  j;  is 
the  system  viscous  dampinq  factor,  kR(y)  is  the  el asto-pl ast ic  spring  restor¬ 
ing  force  function,  k  is  the  system  small  displacement  stiffness,  and  dots 
denote  differentiation  with  respect  to  time.  An  example  af  3  -eaiizat’on  of 


k( y)  is  shown  in  figure  8.  when  the  response  displacement  is  small,  K(y)  is 
linear;  out  when  y  exceeds  the  yield  displacement  U,  the  restoring  force 
becomes  constant,  kD.  This  continues  until  the  sign  of  the  velocity 
changes.  Then  response  starts  to  take  place  along  a  new  straight  line  in  the 
force-displacement  space.  This  continues  until  yielding  occurs  again,  etc. 

As  stated  earlier,  the  response  of  a  linear  SUF  system  to  oroad-band 
random  excitation  is  a  narrow-oand  random  process.  This  means  tnat  wnen 
the  input  has  power  over  a  range  of  frequencies  wider  than  the  bandwidth  of 
the  SDF  system  and  including  its  natural  frequency,  the  system  response  dis¬ 
plays  power  content  mainly  in  the  band  of  frequencies  nearly  surroundina 
The  realizations  of  the  response  stochastic  process  resemble  a  sinu¬ 
soidal  signal  with  randomly  varying  amplitude.  That  is,  the  response  dis¬ 

plays  regularly  occurring  peaks  and  trouqhs  spaced  at  a  time  interval  of 
about  ir/ojg.  This  behavior  also  occurs  in  el asto-pl  ast  ic  systems  when 
yielding  is  not  too  great.  In  fact,  a  half-cycle  'of  response  might  be 
described  b-y  the  curve  shown  in  Figure  9. 

The  starting  amplitude  of  the  response  is  an-_i  and  the  final  amplitude 
is  aj.  Assune  that  a  half-cycle  in  the  response  of  an  el  asto-pl  ast  ic  sys¬ 
tem  is  characteri zed  by  the  curve  shown  in  Figure  9  and  that  this  curve  can  be 
approximated  by 

y(t)  =*  y  (aj_i  +  aj)  ’ t  (aj  ‘  ai_i)  cos  Uf)  t  ,  0  a  <  »/un  (32) 

If  the  dashed  line  in  Figure  9  represents  the  level  at  which  yielding  occurs 
in  the  response,  then  a  portion  of  the  response  during  the  half-cycle  occurs 
in  the  yield  range. 

The  response  can  also  be  depicted  in  the  elasto-plastic  soring  restorina 
force  versus  displacement  space.  When  it  is,  the  system  executes  the  curve 
snown  ’n  -:gure  10. 
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Figure  8.  Force-displacement  curve  for  an  elasto-plastdc  system 
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j-1  Half  Cycle  Number  j 

Figure  9.  Half-cycle  of  narrow-band  response  in  the 
space  of  displacement  versus  time. 


Enerqy  is  dissipated  in  two  ways  during  the  response  half-cycle.  It  is 
dissipated  in  the  elasto-plastic  spring  and  in  the  damper.  The  energy  dissi¬ 
pated  in  the  spring  during  a  half-cycle  equals  the  net  area  under  the  force 
displacement  curve.  This  area  depends  on  the  starting  and  ending  amplitudes, 
and  the  SDF  system  parameters.  The  energy  dissipated  in  the  spring  can  be 
expressed 


7  (aj  ‘  aj-l>  ao  -  D 

7  (°2  *  aj-l}  +  kD(aj  ‘  °)  aj  >  D 


(33) 


This  quantity  clearly  depends  on  the  starting  and  ending  amplitudes  and  on 
the  system  parameters.  The  dissipated  enerqy  can  be  negative,  but  only  when 
the  amplitude  decreases. 

The  energy  dissipated  by  the  damper  during  a  half-cycle  of  response  can 
be  expressed 

ry(’,/«un) 

-  c  /  y  dy  (34) 

M  J  y(0) 


The  response  velocity  can  be  obtained  by  differentiating  y(t)  in  Equation 
32,  and  the  variable  of  integration  can  be  changed  in  Equation  163  to  obtain 


cwnir(  a^ 


aj.l)2/8 


(35) 


The  energy  dissipated  in  the  damper  is  always  positive,  and  depends  on  the 
starting  and  ending  amplitude  of  the  half-cycle  response  and  the  system 
parameters. 
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In  the  analysis  of  Section  II. 2. c  the  half-cycle  transition  probabili¬ 
ties  were  obtained  for  a  linear  system.  In  Section  II. 3. a,  a  method  for 
normalizing  the  transition  probabilities  to  approximately  account  for  all  the 
response  paths  was  developed.  Using  the  methods  established  in  those  sec¬ 
tions,  half-cycle  transition  probabilities  can  be  obtained  for  a' nonlinear 
elastic  system,  based  on  the  analysis  of  Section  II. 4. a.  These  are 


(36) 


Because  of  the  normal ization,  the  sun  of  pA  |A  (ak  |ak  )  over  all  kj 
is  1. 

According  to  the  discussion  given  above,  when  a  nonlinear  system  whose 
force-displacement  diagram  resembles  that  given  in  Figure  10  (and  wnose  damping 
value  is  c)  executes  a  half-cycle  of  response,  it  dissipates  the  energy  given 
by  the  sum  of  Equations  33  ana  35.  The  energy  dissipated  is  a  function  of  a.  . 

J  ^ 

and  aj.’ 

Let  aEj  be  the  random  variable  denoting  the  energy  dissipated  during  a 
half-cycle  of  elasto-plastic,  SOF  system  response.  Then 

=  E0s  +  E0d  <37) 

Since  the  energy  dissipated  in  the  spring  and  damper  depends  on  the  ampli¬ 
tudes  aj_i  and  aj,  and  since  the  amplitudes  of ‘  the  response  are  random, 
AEj  is  also  random.'  The  probability  that  the  realized  value  of  aEj  dur¬ 
ing  a  half-cycle  of  response  falls  within  a  given  interval  is  dependent  upon 
the  joint  distribution  of  Aj_i  and  Aj,  the  random  amplitudes  at  either 
end  of  the  half-cycle.  To  compute  the  probability  distribution  of  energy 
dissipated  as  a  function  of  time,  proceed  in  the  following  manner.  First, 
compute  the  probability  distribution  of  eneray  dissipated  during  each  half- 
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cycle  of  response.  Then  use  this  information  to  compute  the  probability  dis¬ 
tribution  of  energy  dissipated  from  the  start  of  the  response  through  the  jth 
half-cycle  of  response. 

Oefine  a  stochastic  process  Ej,  j  *  0,...,N'  as  the  cunulative  energy 
dissipated  stochastic  process.  It  measures  the  energy  d-issipated  by  a  base 
excited,  elasto-plastic,  SDF  system  at  half  cycles  indexed  zero  through  N'. 
The  cunulative  amount  of  energy  dissipated  through  half  cycle  j  depends  on 
the  amount  of  energy  dissipated  at  half  cycle  j-1  and  the  increment  of  energy 
dissipated  during  the  half- cycle.  This  latter  quantity  is  governed  by  the 
distribution  of  aEj,  j  *  1,...,N'. 

To  compute  the  distribution  of  a£j,  assume  that  during  each  half-cycle 
of  response  the  amplitude  changes  sign.  This  follows  from  the  assumption  of 
narrow  bandedness  of  the  response.  This  means  that  in  the  analysis  only  the 
possibility  that  the  response  will  proceed  from  a  trough  to  a  peak,  or  from  a 
peak  to  a  trough,  is  permitted.  In  Figure  9,  a,  ,  and  a,  will  always  have 

J  *  X  J 

opposite  signs.  The  discrete  amplitudes  defining  possible  values  of  the 
response  are  denoted  aj<,  k  *  1,...,M.  All  the  a^  are  positive  numbers, 
and  when  the  system  goes  through  a  trough  it  simply  assumes  an  amplitude 
whose  value  is  the  negative  of  one  of  the  a^  values.  The  first  of  the 
values  a^,  k  *  1,...,M  where  £  M,  denote  the  response  values  the  system 
can  assume  where  yielding  does  not  occur.  (Recall  that,  for  accuracy,  aM 
must  be  chosen  large  enough  so  that  the  possibility  that  the  response  assumes 
that  value  is  small.)  When  a„,ax  is  the  largest  amplitude  to  be  represented 
in  the  probabilistic  response  analysis,  the  discrete  amplitudes  can  be 
defined  as 


Aa  -  WM 


(37a) 


Using  Equations  33  and  35,  compute  the  maximum  and  minimum  values  the 
dissipated  energy  can  assune  during  a  half- cycle  of  response.  These  are 
denoted  e^ax  and  e^n  and  are  given  by 

'max  *  7  (°2  *  (-^)2)'  k0(4M  -  »)♦  *(*M  -  (-*M1))2/a  ,3*'» 

emin  *  7  (al  "  (*al)2)  +  c“n  w<al  *  (-«!))2/8  (38b) 

where  the  response  at  half-cycle  j-1  Is  required  to  begin  in  the  elastic 
range  and  can  end,  at  half-cycle  j,  as  high  as  a^. 

The  range  of  values  the  dissipated  energy  may  assume  can  now  be  discre¬ 
tized.  The  discretization  increment  is  defined 


ae 


emax  “  emin 
- FT - 


(39a) 


The  discrete  values  that  dissipated  energy  may  assune  can  be  defined 

ek  *  (k  "  7)  Ae  +  em1n  »  k  «  1,...,M'  (39b) 

The  ev  are  uniformly  distributed  in  the  interval  (em^n,  emax). 

Durinq  any  actual  half-cycle  of  response,  the  energy  dissipated  assumes 
one  value  in  a  continuous  range  of  values.  Using  the  present  computation 
scheme  the  range  of  values  is  (emin,  en,ax).  Assume  that  when  the  energy 
dissipated  during  a  half-cycle  of  response  falls  in  the  interval  (ei<  - 
ae/2,  ej<  ♦  ae/2),  it  can  be  accurately  represented  by  e<K.  The  accuracy 
of  this  assunption  is  good  when  M'  is  large  enough. 

Now  recall  that  the  amplitude  transition  probabilities  (Equation  35)  were 
obtained  using  the  same  type  of  discretization  described  above.  The  discrete 
amolitudes  are  a^,  k  3  !,..., M  and  the  amplitude  increment  is  a  a.  When  the 


response  falls  in  the  interval  (aj<  -  aa/2,  a^  +  aa/2),  assume  that  the 
response  is  a^.  For  each  pair  of  amplitudes  -a^,  ko  »  l,...,Mj  at  half- 
cycle  j-1,  and  a^,  ki  *  at  half-cycle  j,  compute  the  range  of  pos¬ 

sible  energies  dissipated  during  the  half-cycle;  the  starting  amplitude  is 
a.  ;  the  ending  amplitude  is  in  the  range  (a.  -  aa/2,  a.  +  aa/2).  The 

K  q  K  2  K  j 

range  of  energy  dissipated  values  is  computed  using  Equations  86  and  88  and 

can  be  denoted  (el"1?  ,  e!"a?  )  where 
<0M  K0K1 


a.  >  D 
*1 


(40a) 


amax 

*k0kl 


-  aa/z)2  -  a^)  -  c.„  .(aki  *  aa/2  -  a*  j2/8  , 
4kJ  *  u,(*k1  *  44/2  '  °)  *  Cun  *(»kl  +  44/2  * 


a.  >  D 
<1 


(40b) 


The  range  of  values  (e™1?  ,  e?13?  )  will  not,  in  general,  eaual  anv  one  of  the 

*0*1 

ranqes  (e^  -  ae/2,  e^  +  ae/2),  k  *  The  probability  pA  |A 

IJ !  j* 1 

a^  ),  of  Equation  39,  can  be  associated  with  the  chance  that  the 

response  will  start  at  an  amplitude  -av  (taken  from  one  of  the  values 

*o 

-a^,  k  *  l,...,Hi)  and  end  at  an  amplitude  a^  ,  with  energy  dissipated  fall¬ 
ing  in  the  interval  (e?1?  ,  e!"a?  ). 

3  k0ki  k<jki 

8ased  on  the  above  arguments,  the  joint,  conditional  pmf  of  amplitude  and 


eneroy  dissipation  increment  can  be  develooed.  Let 


“nn-i'V  kii  V  ’  1  . . 

k2  • 


denote  the  conditional  joint  pmf  of  AEj  and  Aj,  given  Aj_^.  That  is 

P*/j]»J-l(V  \K)  '  ‘  ^2*  AJ  =  ^lIVl  ■ 


kg  * 

k.  *  1, . . . ,M 


k2  *  1, . . . ,M' 


j  *  1 . N  (41) 


The  values  of  the  pmf  are  obtained  in  the  framework  of  a  digital  program  com¬ 


putation  scheme  as  follows.  (1)  All  the  values  of  the  pmf,  p  F  .  i. 

Atj  jlAj>l 

(ek2,  ak j | - ak 0 ) •  are  set  to  zero.  (2)  Set  k0  *  1  and  lq  *  1.  (3)  Compute 

(elc  ok  l*  ek  ok  i  (e“  •  1  and  °AJ|Aj.1(ak1|-ak0>  ■  °AJ|Aj.1(a'|-a')  US,"!I 

Equations  40  and  36.  (4a)  When  the  interval  (eijn,  eTfx)  is  completely 

included  in  one  of  the  intervals  (e^  -  Ae/2,  e^  +  Ae/2),  then  increment  the 
pmf  ^et*  ai  |  ~ai )  by  the  amount  PA  |A  ^a^aj.  (4b)  When  the 

interval  (eu",  e^i  )  overlaps  two  or  more  enerqy  increment  intervals,  say 

(e.  -  Ae/2,  e,  +  Ae/2)  throuqh  (e,  -  Ae/2,  e.  +  Ae/2),  then  the  probabil- 

*1  X1  *2  *2 

ity,  PA  |A  ( a 1 1— a ^)  is  divided  into  fractional  parts  where  the  fractions 

are  linearly  proportional  to  the  fractions. of  the  intervals,  (e,  -  Ae/2, 

*1 

e,  +  Ae/2)  throuqh  (e,  -  Ae/2,  e,  +  Ae/2),  spanned  by  the  interval 

*1  l2  l2 


(e^jn,  e^f*).  These  fractions  of  pA  |A  ^(a^J-a^)  are  used  to  incremeni 
the  values  ^(AjJ  V  "kJ'V  throu9h  »aE,A  |A  /«»,*  Vl'V' 


used  to  increment 


(5)  kj  Is  incremented  by  1,  and  then  return  to  step  (3)  until  kj  =  .  (6)  k, 

Is  Incremented  by  1  and  then  return  to  step  (3)  until  ko  *  M. 

Assune  that  the  starting  probability  distribution  of  response  amplitudes 
at  half-cycle  zero  is  available.  Half  this  distribution  is  denoted  p^Ca^), 
k  *  1,..,,M.  Only  half  the  distribution  is  represented  since  a^  >  0,  k  * 
1, . . . ,M.  The  starting  pmf,  and  all  other  pmf's,  will  be  assumed  symmetric 
about  the  zero  amplitude,  in  what  follows.  In  most  cases  this  pmf  will  be 
zero  for  k  >  Mlt  and  generally 


(0.5  ,  k  -  1 

(0  ,  k  *  1 


(42) 


This  is  the  zero  start  initial  condition.  Two  options  are  available  for 
definition  of  the  starting  pmf,  and  the  formula  given  above  represents  one 
option.  Half  the  probability  can  be  accounted  for,  as  it  is  in  Equation 
42.  In  this  case,  the  marginal  pmf's  and  joint  pmf,  to  be  computed  later  in 
this  analysis,  represent  the  chance  that  specific  amplitude  and  energy 
dissipated  combinations  will  be  realized  during  a  given  half  cycle.  Only  the 
probabilities  for  amplitude  combinations  that  are  negative  at  the  start  of 
the  half  cycle  and  positive  at  the  end  will  be  computed.  8ecause  of  symme¬ 
try,  the  probabilities  of  amplitude  combinations  which  have  positive  ampli¬ 
tudes  at  the  beginning  of  the  half  cycle  and  neqative  amplitudes  at  the  end 
can  be  found  from  those  computed.  Since  the  energy  dissipated  during  a 
response  half  cycle  depends  only  on  the  amplitude  change,  aEj  remains 
unchanged  when  the  signs  on  the  starting  and  ending  amplitudes  are  both 
changed.  (See  Equations  33  and  35.)  The  second  option  for  definition  of  the 
starting  pmf  accounts  for  all  the  probability  by  changing  the  value,  0.5  to 
1.0  in  Equation  42.  In  this  case  A0  must  be  replaced  by  its  absolute  value, 
j  A  o  j ,  because  of  symmetry.  Then,  in  later  comoutations  the  marainal  pmf's 
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'T^1.*-  «  "1M  kj  i ^ 


i 
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and  joint  pmf's  involving  J  Aj  |  are  computed.  As  in  the  previous  case, 
aEj  is  not  affected. 

The  joint  pmf  of  a£j,  the  increment  of  dissipated  energy  during  the  first 
half-cycle,  and  A0  and  Ai,  the  amplitudes  at  the  beginning  and  end  of  the 
first  half-cycle  can  be  obtained  using  Equations  41  and  42.  * 

PaE1A0A1  (Sk2*  V  3,<1^  PaE1AiI  Ao(ek2*  \K 

kg  -  1,  .  .  .  ,Mj 

kj  3  1, . . .  ,M 

k,  *  1, . . . ,M‘  (43) 


Take  PA0^"ak0^  3  pA0^ak0^*  k°  3  l****»Mi  because  of  the  symmetry  assumpti 


ion. 


Only  half  the  complete  joint  pmf  is  represented  here,  but  because  of  the 
synmetry  assumption  the  joint  pmf  evaluated  at  negative  values,  -a^,  is 
the  mirror  image  of  that  evaluated  in  Equation  43.  Specifically, 

kg  —  1, • . * 

p‘WiK'  v  BaEiViK’  "V  \)’  ki  ■ 1 . M 


k~  *  1, . . . ,M‘ 


(43a) 


The  marginal  pmf's  of  aEi  and  Ai  can  be  obtained  by  summing  out  depen¬ 
dence  on  ko  and  klt  then  ko  and  ka,  respectively. 

*  •  p*iK) '  |i  f.i  p*Wi(v  V  \)  ’  k*  * l--"'  (14) 


nl  M 


Pa  (ak  )  *  Z}  2  P*P  a  A  (ek  *  ek  ’  ek  )  ’ 

A1V  V  k^l  k,»l  Ati»W  K2  *0  V 


1 . M 


(*5) 


K-.j 


u 
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When  Equation  45  Is  summed  over  klf  the  result  is  the  value  one-half.  The 
reason  is  that  half  a  symmetric  pmf  is  represented.  When  Equation  44  is 
summed  over  k2,  the  result  is  one-half.  The  reason  is  that,  even  though  all 
the  energy  states  are  represented  (see  Equations  39a  and  b),  only  half  the 
paths  leading  through  these  energy  states  are  considered.  This  fact  is 
accounted  for  by  the  factor  of  1/2  on  the  left-hand  side  of  Equation  44. 

In  the  computation  framework  established  here  it  is  not  feasible  to  allow 
the  system  to  execute  responses  where  yielding  increases  without  limit. 
While  the  transition  probabilities  established  in  Equation  41  can  account  for 
this,  the  computer  storage  required  to  account  for  all  the  paths  would  be 
prohibitive.  Therefore,  a  response  collapsing  procedure  is  established 

here.  This  procedure  concentrates  the  probability  associated  with  all  those 

(m) . 

amplitudes  a^,  k  >  Mlt  into  p^  (a^).  That  is. 


The  (m)  superscript  refers  to  the  fact  that  the  omf  of  At  has  been  modified. 

At  this  point  the  modification  destroys  the  accuracy  of  the  true  ampli¬ 
tude  pmf  at  half-cycle  1.  Indeed.'this  change  affects  all  future  amplitude 
probability  computations.  However,  the  modification  does  not  affect  the 
accuracy  of  the  increment  in  energy  dissipated  probability  computations.  The 
reason  is  this.  Once  a  curve  governing  response  in  the  SOF  system  spr.ing- 
force-versus-alsplacement  diagram  has  oeen  identified,  the  only  quantity 
which  affects  the  enerqy  dissipated  during  a  half-cycle  (besides  the  input 
and  system  parameters)  is  the  starting  amplitude  on  that  curve.  For  example. 
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consider  Figure  11.  If,  at  the  end  of  a  particular  half- cycle  of  response, 
the  system  has  reached  the  yield  threshold,  but  has  not  yielded,  then  durina 
the  next  half-cycle  the  response  will  be  governed  by  curve  No.  1.  If,  on  the 
other  hand,  yielding  occurs  at  the  end  of  a  half-cycle,  then  the  next  half¬ 
cycle  may  be  governed  by  curve  No.  2.  If  the  system  which  yielded  is  excited 
by  the  same  input  as  the  system  wnicn  did  not  yield,  then  the  change  in  dis¬ 
placement  will  be  the  same  in  both  systems,  and  both  systems  will  dissipate  the 
same  amount  of  energy  during  their  response  half-cycles.  This  is  shown  in  in 
tne  figure.  The  conclusion  is  that,  when  a  system  has  yielded  during  a  partic¬ 
ular  half-cycle  of  response,  it  is  accurate,  as  far  as  energy  dissipation 
calculations  are  concerned,  to  treat  that  system  as  though  it  started  its  next 
half-cycle  of  response  from  the  yield  displacement.  This  is  precisely  what  the 
modification  of  Equation  46  does. 

Now  using  Equations  41  and  46  the  joint  pmf  of  aE2,  Al  and  A2  can  be 
computed. 


k2  *  1, 


.M' 


(471 


This  can  be  used  to  obtain  the  marginal  pmf's  of  4E2  and  A2.  Then  the  mar¬ 
ginal  pmf  of  A2  can  be  modified  as  in  Equation  46.  Then  the  joint  pmf  of 

&E3,  A2  and  A3  can  be  obtained,  and  so  on,  until  all  the  pmf's  paE  (e^), 

J 

J  *  1,...,N' ,  k  »  1, . . . ,M'  are  known. 

Finally,  the  specification  for  the  stochastic  process  E j ,  j  =  0,...,N‘, 

can  be  formed  using  the  pmf's'  p^£  (ek),  j  *  1,...,N',  k  =  1,...,M'.  Ej  is 

j 

the  energy  dissipated  by  the  system  throuan  half-cycle  j 


n  'j' 

m  y 


Displacement 


Force  displacement  response  of  elasto-plasti 
SDF  systems. 


Is  a  Markov  process.  Hie  energy  dissipated  by  the  system  through  naif-cycle  j 
depends  only  on  the  energy  dissipated  through  half-cycle  j-1  and  the  energy 
dissipated  Increment  in  the  half-cycle  between  j-1  and  j. 

Let  p-  ic  (e.  |e.  ),  j  s  1,...,N',  k0,  k x  3  1,...M",  be  the  conditional 

j  •  j  - 1  kl1  k° 

pmf  for  the  random  process,  Ej,  j  *  0....N',  defining  transition  probabili¬ 
ties  between  states  of  dissipated  energy.  Let 

M"  *  N'  x  M*  (48) 


ek  »  (k  -  f)  ae  ,  k  a  1, • • • ,M" 


where  ae  is  defined  in  Equation  39a. 

Then  the  transition  probabilities  are  defined 


PEjiEj-iKieko)'PiEjK'tko) 


j  =  1 . N' 


That  is,  the  chance  that  the  dissipated  energy  state  changes  from  e.  to  e. 

Ko  Ki 

in  one  half-cycle  equals  the  chance  that  eu  -  e.  units  of  energy  are  dissi- 

K  i  Kg 

pated  during  the  jth  half-cycle. 

Let  pc  (ev),  j  *  0,...,N',  k  =  1, . . . ,MM,  be  the  pmf  of  enerqy  dissipated 

j 

through  the  jth  half-cycle.  Then  the  joint  pmf  of  energy  dissipated  throuah 
the  j-lst  and  jth  half-cycles  is  obtained  following  Equation  34a;  it  is  ■ 


Pej-iej(V  \)  ‘  Vj-iKl'O  Pej-iK)’ 


j  -  1 . N-, 

kQ,k1  »  1, . . . ,M" 


The  marginal  pmf  of  energy  dissipated  through  the  jth.  half-cycle  is 


Pr  (c k  )  a  23  Pr  f  (ck  *  Ek  ) 

Ej'V  kn»l  Ej-lEj\  k0  V 


j  3  1,...,N’ 
k  .  -  1 . M" 


>vv>v->y- v-y-y-is.  vS*. 


Or  by  combining  the  operations  In  Equations  51  and  52,  the  marginal  pmf  of  E. 
from  the  pmf  of  E^  can  be  obtained 


.i  3 

(53) 


By  specifying  the  starting  pmf  of  Eo;  one  can  obtain  the  pmf  of  Ej  for 
any  j  *  through  successive  applications  of  Equation  53. 

Note  that  no  assumptions  regarding  the  stationarity  of  the  input  signals 
has  been  made  in  the  analysis  of  this  section;  therefore,  the  probability 
distribution  of  the  enerqy  dissipated  response  excited  by  a  nonstationary 
input  can  be  computed  using  the  technique  developed  above.  Further,  energy 
dissipation  caused  by  Inelasticity  in  the  spring  and  viscous  damping  need  not 
be  considered  simultaneously.  No  feature  in  the  theoretical  development  of 
this  section  precludes  the  consideration  of  energy  dissipation  due  to  one 
source  only. 

When  the  pmf,  pF  (e.  ),  is  summed  over  all  kj  the  result  will  always  be  a 

j 

value  of  1.  The  reason  is  that  all  potential  values  the  dissipated  energy 

might  assume  have  been  approximately  accounted  for.  The  definition  of  e,,,^ 

and  e  in  Equations  38a  and  38b,  and  use  of  these  in  Equation  39a  and  39b, 
min  ’ 

the  definition  of  the  dissipated  energy  states,  has  guaranteed  this. 

The  accuracy  of  computations  performed  using  the  approach  specified  above 
will  improve  as  the  values  of  M,  Mlf  M'  and  are  increased.  The  first 
three  values  involve  the  fineness  of  the  discretization;  the  last  value  is 
related  to  the  displacement  interval  where  response  paths  are  accounted  for. 


b.  Numerical  example--This  section  gives  a  numerical  example  demonstrating 


use  of  a  computer  program  which  implements  the  analyses  presented  in  Section 


—  tTt-'  ■.-»  iv1".  WA  "TTT  u  '  V*  WW.WWVTV 


II. 5. a.  In  this  example  the  probability  distribution  of  energy  dissipated 
up  to  a  particular  time  is  computed  as  a  function  of  time.  Specifically,  the 

prof.  Pc  (O.  J  *  0,...,N't  k  *  1,...M",  is  computed  at  the  half-cycle  points 

\i  * 

of  an  SOF  system  response.  The  input  used  to  excite  the  system  is  a  band- 
limited  white  noise.  The  system  computation  and  input  parameters  are  listed 
in  Table  4. 


TABLE  4.  SYSTEM,  COMPUTATION,  AND  INPUT  PARAMETERS  FOR 
ENERGY  DISSIPATED  pmf  CALCULATION 


At  *  ir/CSu^) 

o'x'  =  1.0 

C  =  0.05 
M  *  18 
Mx  *  12 
M'  *  20 
M"  *  1000 
a,^  *  15.0 

N'  *  50 


=  6.28 


k  »  39.5 


0  *  10.0 


The  response  pmf  was  computed  at  a  time  interval  of  ten  half-cycles.  The 
resulting  pmf’s  were  interpolated  so  that  the  results  could  be  displayed  as  a 
sequence  of  pdf's.  The  results  are  shown  in  Figure  12.  Each  pdf  is  plotted 
in  a  vertical  plane  of  pdf  ordinate  versus  energy  dissipated. 

At  time  zero  the  system  starts  with  no  energy  dissipated,  so  the  first 
pdf  is  a  delta  function.  As  time  progresses,  the  chances  for  increased 
amounts  of  dissipated  energy  increase,  and  the  pdf's  of  Ej  spread.  The 
pdf's  shown  In  Figure  12  resemble  exponential  Ddf’s,  and  simple  computations 


show  that  the  ’-esults  are  nearly  exponentially  distributed. 
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The  mean  and  variance  of  each  pdf  were  evaluated  by  the  computer  pro¬ 
gram.  The  graph  of  the  mean  value  of  energy  dissipated  versus  time  is  plot¬ 
ted  in  the  horizontal  plane  of  Figure  12.  This  curve  increases  with  time  and 
appears  to  be  in  a  condition  of  constant  (straight  line)  increase  past  half¬ 
cycle  30. 

A  listing  of  the  computer  program  used  to  perform  this  numerical  example 
is  given  in  the  Appendix. 


6.  FIRST  PASSAGE  OF  DISSIPATED  ENERGY  IN  AN  ELASTO-PLASTIC  SYSTEM 

a.  Theoretical  analysis — The  analysis  of  Section  II. 5. a  pointed  out 

that  the  sum  of  the  marginal  pm f  of  dissipated  energy,  pc  (c.  ),  j  «  0, ....N1, 
k  *  1 , . . . ,MM ,  over  all  k  equals  1,  because  all  the  values  that  dissipated 
energy  can  assume  were  accounted  for.  Also  pointed  out  was  tnat  the  sto¬ 
chastic  process  tracking  dissipated  energy  is  Markov.  The  first  passage 
problem  for  energy  dissipated  in  an  elasto-plastic  system  can  be  solved  by 
taking  advantages  of  the  Markov  property  and  rearranging  the  computation 
scheme  to  permit  the  escape  of  some  response  paths  past  a  barrier. 

Since  the  stochastic  process  tracking  dissipated  energy  is  Markov,  its 
transition  probabilities  can  be  arranged  in  a  transition  matrix  defined  as 
f o 1 1 ows . 


Cjnl  * 


PE j+1 1  Ej(ell  eM" )  PE j+ L |  E j ( e2  lxM" ) 


PE  j+1 1 E  j  ^ e MM  I  eM"  ^ 


Each  of  the  elements  in  the.  matrix  is  defined  as  in  Equation  50.  And  when 
the  pmf  of  dissipated  energy  is  arranged  in  vector  form, 

{pj}  *  (pE  (ej)  PE  (e2)  •••  PE.(eMH))T  ,  j  *  0,...,N'  (55) 

J  J  J 

Through  specification  of  the  pmf  of  dissipated  energy  at  half-cycle  zero, 
denoted  {p0},  one  can  obtain  the  pmf  of  dissipated  energy  at  half-cycle  j. 
This  is 

i  ■  A  (p0)  <m> 

(Details  of  this  general  procedure  are  given  in  Section  II.2.d.) 

The  computation  procedure  outlined  above  can  be  modified  to  yield  first 

★ 

passage  probabilities.  Let  [Pj+^J  denote  a  square  submatrix  of  [Pj+j]  taken 
from  the  upper  left  corner.  Then 


_  ★  - 

fj.ii  ■ 


j  •  (57) 

where  L  <  M".  And  let  {po}  denote  a  vector  of  elements  to  be  operated  on  by 
the  Cpj}.  The  elements  in  {po}  are  the  probabilities  that  Eo  occupies  the 
states  ei  through  e^. 

(PnJ  *  f°r  (eO  pc  Uo)  ***  Pc  (e.))T 


(58) 


The  star  (*)  superscript  is  attached  as  an  indicator  that  the  transition  and 
state  probabilities  corresponding  to  the  entire  range  of  realizations  are  not 
included  in  Cpj+j]  and  {p 0} -  A  vector  { p^ }  can  be  obtained  through  succes¬ 
sive  operations  of  [p£]  on  {po}.  This  is 


fell  "  px 


j  -  1 . r 


(59) 


f  ★  - 

The  elements  of  {p. }  represent  the  probabilities  that  particular  states  in 

j 

the  energy  dissipated  range,  e;  through  «{_•  be  reached  following  paths 

that  do  not  pass  outside  the  level  e[_.  Therefore,  at  half-cycle  j,  the  no 
passage  probability  of  dissipated  energy  is 


P(Tj  >  j*/«n) 


PE 


j  *  1 . N' 


(60) 


T x  is  the  random  variable  denoting  the  time  at  which  first  passage  of  the 
dissipated  energy  response  passes  outside  the  level  ej_.  »/u»n  is  the 

duration  of  a  half-period  of  the  system  under  consideration,  p*  (e^)  is  the 

* .  J 

kth  element  in  the  vector  {p^}. 

The  first  passage  event  is  the  complement  of  the  no  passaae  event;  there¬ 
fore,  the  first  passage  probability  is  given  by 

P(T1  <.J*/^)  *  1  -  P(TX  >  Jt/mn)  ,  j  «  1 . N*  (61) 


This  first  passage  probability  is  a  nondecreas inq  function  of  j. 

As  long  as  the  input  excitinq  response  in  Equation  31  has  nonzero  mean 

square  power,  the  first  passage  probability  will  increase  as  a  function  of 
time.  This  corresponds  to  a  feature  in  the  matrix  of  Equation  57  where  the 
rows  ao  not  aaa  to  1;  that  is,  the  probability  that  the  response  starts  at 
any  point  in  the  range  through  ^  and  passes  to  another  point  in  the  same 
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range  Is  lower  than  1.  The  remainder  of  the  probability  corresponds  to  paths 
which  pass  outside  the  barrier  e^. 

The  analysis  done  in  this  section  is  restricted  by  the  assumptions  used 
in  obtaining  the  elements  in  the  transition  matrix.  Equation  54.  Most 
Important,  the  response  is  assumed  narrow  band;  therefore,  the  analysis  is 
accurate  only  for  systems  In  which  large  amounts  of  yielding  cannot  occur 
during  most  half-cycles.  First  passage  of  energy  dissipated  due  to  the  non¬ 
stationary  response  of  an  SDF  system  can  be  accounted  for,  though,  since  the 
probability  distribution  of  energy  dissipated  in  nonstationary  response  can 
be  computed  using  the  underlying  approach  of  Section  II. 5. a. 

b.  Numerical  examples — This  section  sunmarizes  some  numerical  examples 
demonstrating  the  use  of  a  computer  program  which  implements  the  analysis  of 
Section  II. 6. a.  The  computer  program  used  here  involves  only  slight  modifi¬ 
cations  in  the  computer  program  of  Section  II. 6. b;  therefore,  it  is  not 
listed  in  the  Appendix.  The  first  passage  probability  functions  are  computed 
for  three  cases.  In  each  case  the  input  excitation  is  a  band-limited  white 
noise  stochastic  process.  Three  barrier  levels  are  chosen  for  analysis.  The 
system,  computation  and  input  parameters  are  listed  in  Table  5. 

The  first  passage  probability  was  computed  at  an  interval  of  every  ten 
half-cycles.  The  first  passage  probability  results  were  interpolated  so  that 
they  could  be  displayed  as  continuous  curves.  These  are  shown  in  Figure  13. 

A  zero  start  condition  was  used  on  each  first  passage  computation;  that  is, 
the  system  was  started  with  zero  energy  dissipated.  For  this  reason,  each 
curve  starts  at  zero.  As  time  progresses  the  chance  that  the  response  passes 
outside  the  barrier  increases,  so  the  first  passage  curve  increases.  As  time 
goes  to  infinity  these  curves  approach  one  since  the  input  is  stationary; 
excursion  of  any  finite  barrier  is  assured  as  time  goes  to  infinity.  As  the 
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TABLE  5.  SYSTEM,  COMPUTATION,  AND  INPUT  PARAMETERS  FOR  FIRST  PASSAGE 
OF  DISSIPATED  ENERGY  CALCULATION 
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barrier  height  is  increased  the  first  passage  probability  tends  to  increase 
more  slowly  with  time. 


7.  PROBABILITY  OF  FAILURE 

a.  Theoretical  analysis— The  reliability  of  a  structure  is  the  probabil¬ 
ity  that  it  will  perform  satisfactorily  over  a  preestablished  period  of 
time.  The  requirement  for  satisfactory  performance  implies  that  some  crite¬ 
rion  has  been  established  to  judge  whether  or  not  the  response  is  accept¬ 
able.  Failure  of  a  system  to  satisfy  the  response  criterion  is  called  struc¬ 
tural  failure  and  the  probability  of  failure  is  1  minus  the  reliability. 

In  many  situations  failure  occurs  only  with  the  physical  collapse  of  a  struc¬ 
ture;  and  in  such  cases  the  reliability  is  the  chance  that  no  collapse  will 
occur  during  the  design  life  of  a  structure.  When  structural  collapse  is  the 
failure  criterion  the  reliability  is  difficult  to  estimate  accurately  since  a 
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structural  analysis  predicting  collapse  must  be  executed  in  the  reliability 
analysis.  For  this  reason  structural  reliability  is  often  simply  bounded  or 
estimated  approximately,  using  approximate  nonlinear  structural  analyses 

where  necessary. 

There  are  three  features  which  may  lead  to  randomness  in  the  dynamic 
response  of  a  structure.  These  are  randomness  in  (1)  the  input,  (2)  the 

mechanical  properties  of  the  structural  material,  and  (3)  the  geometry  of  the 
system  under  consideration.  When  the  material  properties  and  geometry  of  a 
structure  are  practically  deterministic,  then  only  the  input  is  random;  and 
the  reliability  analysis  can  be  reduced,  in  many  cases,  to  a  first  passaqe 

analysis  like  the  one  presented  in  Section  II. 2  or  II. 3  or  II. 6.  The 

deterministic  failure  level  of  the  system  is  simply  chosen  as  the  first  pas¬ 
sage  barrier,  and  the  failure  probability  is  computed  as  the  first  passage 
probability. 

In  a  more  general  situation  the  mechanical  properties  of  the  structural 
material  are  random  and  no  single  first  passage  analysis  can  be  used  to 
estimate  the  failure  probability.  Rather,  a  combination  of  information 
regarding  the  probabilistic  character  of  the  material  must  be  used  with 
information  on  the  probabilistic  character  of  the  peak  response  to  estimate 
the  chance  of  failure.  Specifically,  through  testing,  the  conditional  proba¬ 
bility  of  failure  of  a  material  can  be  developed.  The  probability  of  failure 
is  conditional  on  some  measure  of  the  load  on  the  material.  For  example, 
when  a  test  specimen  is  loaded  monotonical ly  the  probability  of  failure  is 
conditioned  on  peak  stress  or  peak  strain.  When  a  test  specimen  is  loaded 
cyclically  the  probability  of  failure  may  be  conditioned  on  a  measure  of 
accumulated  plastic  deformation,  or  energy  dissipated,  etc.  The  following 
material  considers  this  problem  for  SOF  systems.  The  response  of  an  elasto- 
plastic  system  is  considered,  and  the  system  under  consideration  is  assumes 


to  have  a  random  capacity  for  dissipating  energy.  Through  testing,  the  prob¬ 
ability  of  failure  can  be  specified,  conditional  upon  the  amount  of  enerqy 
dissipated  in  the  system.  When  the  letter  F  is  used  to  denote  the  failure 
event,  the  conditional  probability  of  failure  is  denoted  P(f|ec  ■  e),  e  £ 
0.  denotes  the  total  energy  dissipated  over  the  entire  duration  of  the 
system  response,  and  the  assumption  is  made  that  the  conditional  probability 
of  failure  is  known  for  all  realization  of  dissipated  energy,  e.  When  the 
probability  distribution  of  E  is  specified,  then  this  can  be  used  with 

w 

P(F|Ec  «  e)  to  find  the  probability  of  failure. 

The  quantity  Ect  whose  probabilistic  character  needs  to  be  specified,  is 
the  peak  value  of  energy  dissipated  over  the  entire  response  duration.  The 
cdf  of  this  quantity  can  be  found  in  terms  of  the  first  passage  probability 
for  dissipated  energy  obtained  in  Section  II. 6.  The  following  material 
slightly  changes  the  notation  established  in  Equations  60  and  61.  Let  T 

c 

denote  the  random  variable  representing  the  time  at  which  first  passage  of 
the  dissipated  energy  barrier  level  e  occurs.  If  Tc  is  the  duration  of  the 
response,  then  P(TC  £  Tc)  is  the  probability  that  first  passage  beyond 
the  barrier  e  occurs  at  or  before  Tc,  and  can  be  found  from  Equation  61. 

But  this  probability  can  be  viewed  in  another  way.  It  is  the  chance  that 
the  peak  value  of  energy  dissipated  in  the  time  interval  (0,  Tc)  surpasses 
the  value  e.  That  is 

P(EC  >  e)  =  P(Te  £  TJ  (62) 

The  complement  of  the  event  whose  probability  is  described  on  the  left-hand 
side  is  that  the  energy  dissipated  is  less  than  or  equal  to  e;  therefore, 

P(EC  <  e)  «  1  -  P(Te  £  Tc)  (S3) 


The  function  on  the  left  can  be  computed  through  repeated  application  of  the 
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first  passage  analysis.  That  is,  a  first  passage  analysis  can  be  run  through 
time  T  ,  for  various  values  of  the  barrier  level,  e;  in  this  way  P(EC  £  e) 
is  evaluated  at  a  discrete  set  of  values. 

The  pdf  of  Eg  can  be  obtained  by  differentiating  the  cdf  of  E,-,  Equa¬ 
tion  63,  with  respect  to  e.  The  result  is 

fE  (e)  s^P(Ec<0  ,  e  >  0  (64) 

The  product,  fP  (e)  de,  denotes  the  probability  that  E  has  a  realiza- 
fcc 

tion  in  the  interval  (e,  e  +  de];  therefore,  the  joint  probability  that  sys¬ 
tem  failure  occurs  and  the  dissipated  energy  lies  in  the  interval  (e,  e  +  de] 
is  given  by 

P(F,  e  <  Ec  <  e  +  de)  -  P(F|EC  *  e)  f£  (e)  de  ,  e  >  0  (65) 

The  probability  of  failure  can  be  obtained  by  integrating  out  dependence  on  e 
in  the  above  expression.  The  result  is 

P(F)  P(F|EC  -  e)  f^e)  de  (66) 

The  above  expression  can  be  discretized  so  that  the  discrete  nature  of  the 
cdf  of  Ec  can  be  directly  accounted  for. 

The  computations  performed  above  assume  that  the  SOF  system  will  fail  due 
to  the  dissipation  of  energy  in  the  system.  In  fact,  failure  may  be  more 
strongly  related  to  another  measure  of  structural  response  or  a  collection  of 
other  measures.  If  the  conditional  probability  of  failure  given  these  other 
measures  of  response  can  be  specified,  either  based  on  experiments  or  theore¬ 
tical  considerations,  and  if  the  probabi 1  istic  character  of  these  other  mea¬ 
sures  of  structural  response  can  be  obtained,  based  on  the  type  of  analysis 
given  in  this  report,  then  the  probability  of  failure  can  be  specified  using 
the  same  general  approach  developed  in  this  section.  In  view  of  this,  the 


failure  analysis  presented  in  this  section  is  simply  meant  to  demonstrate  a 
technique  which  takes  into  account  the  randomness  in  the  excitation  and  the 
system  characteristics. 


b.  Numerical  example — This  section  presents  a  nunerical  example  demon¬ 
strating  the  use  of  a  computer  program  which  implements  the  failure  analysis 
of  Section  II. 7. a.  The  computer  program  uses  the  first  passage  analysis  of 
Section  II. 6. b  to  develop  the  peak  response  cdf;  therefore,  the  input  used 
in  that  section  is  in  effect  here.  That  is,  a  white  noise  is  used  to  excite 
the  system.  The  time  duration  of  the  input  and  response  are  taken  to  be  50 
half-cycles  of  the  response.  The  cdf  and  pdf  of  ^  are  evaluated  at  that 
time.  The  system,  computation,  and  input  parameter  used  to  obtain  the  cdf  of 
Ec  are  given  in  Table  6. 

TABLE  6.  SYSTEM,  COMPUTATION,  AND  INPUT  PARAMETERS 
FOR  PROBABILITY  OF  FAILURE  CALCULATION 

At  »  ir/(5un) 

*  1.0 

C  =*  0.05  k  »  39.5  0  *  10.0 

M  *  18  Mj  ■  12  M*  a  20 

aM  3  15.0 

N'  a  50 

cdf  of  Ec  evaluated  at 
e3  *  840  AE  »  280 

c4  *  1120 
e5  »  1400 

Eg  *  1680 

e7  «  1960 
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Part  of  the  cdf  of  Ec  was  computed  at  increments  graphed  in  Figure  14. 
The  cdf  was  computed  at  increments  of  dissipated  energy  equal  to  280.  The 
cdf  is  shown  as  a  step  function.  The  pdf  is  estimated  by  differencing  the 
cdf.  The  increments  obtained  by  differencing  the  cdf  are  shown  in  Figure  15. 

In  order  to  compute  the  failure  probability,  the  conditional  probability 
of  failure  given  the  level  of  dissipated  energy  must  be  known.  In  this  exam¬ 
ple  it  is  assumed  that  the  conditional  probability  of  failure  is  that  shown 
in  the  graph  of  Figure  18.  The  dissipated  energy  increments  upon  which  the 
failure  probability  is  conditioned  were  chosen  to  correspond  to  those  used  in 
obtaining  the  cdf  of  Ej.  This  was  done  for  convenience,  though  any  discre¬ 
tization  in  the  conditional  failure  probability  function  could  be  interpo¬ 
lated  to  make  it  compatible  with  the  increments  in  the  cdf  of  Ec. 

The  discretized  form  of  Equation  66  used  to  obtain  the  failure  probabil¬ 
ity  is 


P(F) 


•  x  £  p<pi 


i 


Ec  '  ?(ek  *  W>  tp<Ec  <  W  -  p<Ec 


where  AE  is  the  dissipated  energy  increment  e^,  k  =  1, ...,M2 
pated  energy  levels  where  the  cdf  of  Ec  is  computed,  and  M2 
of  dissipated  energy  values  where  the  computation  is  carried 
oresent  example,  the  probability  of  failure  assumes  the  value 


<  ek)]  (67) 

are  the  dissi- 
is  the  number 
out.  In  the 


P(F)  *  0.0(0.6057)  +  0.2(0.1979)  +  0.5(0.1031)  +  0.8(0.0532)  +  1.0(0.0401) 

=  0.1750  (68) 

This  example  shows  that  it  is  necessary  to  obtain  the  cdf  of  Ec  over 
these  values  only  where  the  conditional  probability  of  failure  is  greater 
than  zero  ana  less  than  1.  None  of  the  energy  dissipated  quantities  can  lead 
to  failure  when  the  conditional  probability  of  failure  is  zero,  and  all  tne 

energy  dissipatea  quantities  lead  to  failure  wnen  the  conditional  probability 
of  failure  is  1. 
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Figure  15.  Approximate  partial  pdf  of  energy  dissi 
pated  during  50  half  cycles  of  response 
of  an  elasto-plastic  SDF  system. 
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Figure  16.  Conditional  probability  of  failure  of  an 
elasto-plastic  SDF  system  given  dissipated 
energy. 


In  summary,  this  section  has  shown  an  approach  to  the  prediction  of  reli¬ 
ability  of  SUF  systems  that  are  nonlinear  and  subjected  to  both  stationary  and 
nonstationary  loads.  Various  measures  of  response  have  been  considered  and  the 
failure  probability  has  been  computed  for  a  system  with  a  random  failure  level 


of  dissipated  energy. 


III.  FAILURE  ANALYSIS  OF  MULTIPLE-DEGREE-OF-FREEDOM  SYSTEMS 


1.  INTRODUCTION 

In  the  previous  section,  a  probabilistic  analysis  of  single-degree-of- 
freedom  (SUF)  systems  was  presented.  The  extension  of  the  problem  to 
multiple-degree-of -freedom  (MOF)  systems  is  not  trivial.  Recall  that  one  of 
tne  disadvantages  of  the  previously  presented  SOF  approach  was  the  storage 
requirements  for  the  state  probabilities  and  the  transition  probabilities. 

Even  when  some  memory  requirements  were  reduced  by  additional  computational 
effort,  tne  memory  storage  was  large.  Any  solution  for  MDF  systems  should 
include  considerations  for  a  reasonably  large  number  of  degrees  of  freedom  (at 
least  several  hundred).  This  constraint  has  led  to  an  alternate  procedure  for 
predicting  the  probability  of  survival  of  structures  excited  by  highly  tran¬ 
sient  loads,  such  as  blast  and  shock. 

The  survival  probability  of  a  structural  system  subjected  to  blast  and 
shock  loads  can  be  computed  using  the  basic  elements  described  in  the  opening 
section  of  this  report.  Tne  random  character  of  the  structural  system  and  of 
the  loadings  must  be  known.  A  mathematical  procedure  must  be  developed  that 
models  the  system  so  tnat  the  random  response  is  predicted.  The  failure  level 
must  be  char eicteri zed  and,  finally,  intersection  of  response  and  failure  must 
be  integrated  into  a  probability  of  failure  or  its  complement,  the  probability 
of  survival.  This  section  presents  an  approach  using  these  elements  for 
Inelastic  systems  that  are  represented  by  MDF  solutions  that  have  some  chance 
of  failing  at  any  of  a  number  of  locations. 

2.  THEORETICAL  ANALYSIS 

a.  Failure  at  a  single  point--The  survivable  probability  of  a  blast- 
excited  structure  can  be  computed  by  a  four-step  procedure:  (I)  A  numerical 
scneme  for  computing  the  response  at  various  locations  on  an  inelastic 


structure  is  selected.  (2)  The  parameters  of  the  Input  and  the  structure  are 
characterized  probabilistically.  (3)  Peak  responses  are  expanded  in  a  series 
involving  the  input  and  system  parameters,  and  the  mean  and  variance  of  the 
peak  response  at  each  location  are  determined.  (4)  Margins  of  survival  are 
established  and  used  to  compute  the  probability  of  survival.  Each  of  these 
steps  will  now  be  described  in  detail. 

The  equation  governing  the  response  of  an  MDF  structure  can  be  written 

[m]  {z}  +  { R( z)}  *  { F( t ) }  (69) 

where  [m]  is  the  mass  matrix,  {z|  is  the  displacement  vector,  dots  indicate 
differentiation  with  respect  to  time,  (R(z)}  is  the  restoring  force  vector, 
ana  {F(t)|  is  the  forcing  vector.  When  the  system  characteristics  are  random, 
the  parameters  governing  the  behavior  of  {R(z)}  are  random  variables.  When 
the  input  is  a  random  process,  the  parameters  of  tne  forcing  vector  are  random 
variaoles. 

Let  8^.,  j=l,  •••,  m,  denote  the  random  structural  parameters.  Let  8^, 
j=m+l,  •••,  n,  denote  the  random  input  parameters,  then  the  response  at  a 
point,  i,  on  the  structure  is  a  randan  process  and  can  be  expressed 

*  z.  (t,  gj,  •••,  8n).  A  number  of  techniques  can  be  used  to  compute  this 
response  when  the  3j,  j*l,  •••,.n,  are  specified.  In  reality  most  dynamic 
problems  do  not  explicitly  solve  the  system  of  equations  shown  in  Equation  69, 
but  rather  the  system  is  spatially  discretized  in  some  way  and  the  equations 
of  motion  are  solved  on  a  cell-by-cell  basis  In  a  sweep  through  the  system  for 
a  small  time  step.  Simplified  approaches  can  be  used  as  well  as  complex 
finite  element  solutions. 


(70) 


t  h 

The  maximum  response  in  time  at  the  i  point  is 


zi(*i»  5n)  »  mJx  z^t,  8lt  •••,  Sn) 


Note  that  although  Zi  and  z^  represent  displacement  here,  they  could  represent 
any  measure  of  the  response  and  its  maximum.  Let  ik  and  o?  represent  the  mean 
and  variance  of  the  random  variable  8j,  j»l,  ••*,  n,  and  let  p<t  be  the  corre¬ 


ct 

lation  coefficient  for  8j  and  8^.  The  function  Z  can  be  expanded  in  a  Taylor 
series  about  the  means  of  the  parameters. 


n  3Zi 

[»l  *  Pi  "i 


l»l 


(8j  -  Uj)  + 


(71) 


where  {u}  is  the  vector  of  mean  values  of  the  8^,  j«l,  •••,  n. 

The  mean  and  variance  of  the  peak  response  can  be  computed  approximately 
as 


Wt'lM 


n  n  3Z , 

V[Zt]  ^  £  £  — L 

j-1  k-1  38j 


3Zj 


{»} 
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{u} 


pjk  °j  \ 


(72a) 


(72b) 


These  expressions  are  approximate  because  higher  order  terms  in  the  series 
expansion  have  been  neglected.  These  terms  can  be  retained  when  necessary. 
Several  appropriate  evaluations  of  the  importance  of  the  higher  order  terms 
can  be  considered.  If  higher  order  terms  do  not  significantly  change  the 
solutions,  they  need  not  be  retained.  Also  if  the  peak  response  function  can 
be  easily  differentiated,  either  analytically  or  nunerical ly,  the  slope, 
curvature  or  higher  order  differentials  can  be  evaluated  on  either  side  of  the 
mean  to  detect  changes.  If  the  changes  are  zero  or  small,  no  further  terms 
should  be  necessary.  An  equivalent  approach  is  to  examine  the  order  of  the 


peax  response  function  in  a  reasonable  range  of  the  mean  determined  by  the 
failure  level. 

Structural  failure  is  assuned  to  occur  when  tne  peak  response  at  any 
single  point  exceeds  the  failure  level.  Let  be  the  random  variable  denot¬ 
ing  tne  failure  level  at  a  point  i.  Let  E[L^]  and  V[L.j]  be  the  mean  and  vari¬ 
ance  of  L^.  Assume  that  and  L^  are  independent.  A  margin  of  survival  at 
point  i,  M . ,  can  then  be  established  as 

Mi  -  Li  -  Zi  (73) 

This  margin  has  mean 

EOi]  -  E[Li]  -  E[Zi]  (74) 

ana  variance 

V[Mi]  -  V[Li]  +  V[Zi]  (75) 

* 

Failure  occurs  when  the  margin  of  survival  is  negative  and,  therefore,  the 
probability  of  survival  is 

Fs  ■  P(rti  >  0)  (76) 

Tnis  can  usually  oe  evaluated  using  E[M.]  and  V[M^]  when  tne  distribution 
of  Mi  is  specified.  When  Zi  and  Li  are  either  normally  or  lognormally  dis- 
trioutea,  Ps  can  easily  be  determined. 

A  development  analogous  to  the  above  can  be  executed  when  Z  is  any  mea¬ 
sure  of  tne  structural  response. 

The  following  numerical  example  considers  the  case  where  tne  potential 
for  failure  exists  at  only  one  point  on  a  structure. 
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b.  Numerical  example— A  laterally  supported  vertical  cylinder  composed 
of  two  materials  can  oe  simply  modeled  as  a  set  of  axial  springs.  The  two 
materials  are  modeled  as  springs  in  parallel,  while  the  structure  length  is 
modeled  as  springs  in  series,  An  example  structure  and  model  are  shown  in 
Figure  17.  The  materials  are  assumed  elasto-plastic. 

Consider  the  structure  excited  by  a  random  blast  input  defined  by 
P  »  a  +  ae"at  where  P  is  the  pressure  history.  The  peak  pressure  A'  *  A+B  is 
a  normally  distributed  random  variable  with  mean  1280  psi  (8.8  MPa)  and  stan¬ 
dard  deviation  102  psi  (0.7  MPa),  8  *  -1.069A,  and  o  is  a  normally  distributed 
random  variable  with  mean  0.80  and  standard  deviation  0.2  t  as  time. 

The  peak  response,  taken  here  as  strain,  at  point  C  (Fig.  17)  has  been 
predicted  for  the  mean  values  of  the  input  using  a  numerical  approximation  to 
the  solution  of  Equation  69  (Ref.  53).  Additionally  the  peak  responses  near 
the  means  provide  partial  differential  quotients  to  approximate  the  partial 
differentials  in  Equation  72b.  The  results  obtained  are 

z|{uj  *  2.12  (.10-3); 

3Z/3A|^  »  0.429  (lO-3); 

3Z/3a  =•  -0.300  (lO”3); 

and  p  ,  is  assumed  0.1. 

art 

Using  Equations  72a  and  b,  the  mean  and  variance  of  Z  become  2.12  ( 10“3 ) 
and  9.0  ( 10-8) .  Let  L  be  the  strain  at  which  failure  occurs  and  define  the 
mean  of  L  as  0.002  and  2.5  (10-7),  respectively.  From  the  above  development, 
tne  mean  and  variance  of  M,  the  margin  of  survival,  become  um  =  1.2  (10-1*)  and 
V  »  3.4  (l0-;). 


Assuming  a  normal  distribution  for  Z  and  L  gives  the  probability  of  sur 
vival  as 


Ps 


1  -  + 


1.2  (10-**) 
5.8  (10-4) 


41.8% 


where  is  the  standard  normal  cumulative  distribution  function. 

In  the  present  application,  sufficient  accuracy  was  obtained  using  the 
first  two  terms  in  tne  Taylor  series  (Equation  71).  More  complicated  situa¬ 
tions  may  require  inclusion  of  additional  terms. 

c.  Failure  at  multiple  polnts~~Note  the  possibility  of  multiple  point 
failure.  That  is,  the  probability  of  structural  failure,  given  that  it  could 
conceivably  fail  at  any  of  a  number  of  locations,  N.  If  the  failure  level  at 
one  location  is  independent  of  failure  at  other  locations,  then  the  solution 
is  simple.  If,  however,  the  failure  level  at  one  location  is  correlated  to 
failure  at  other  locations,  the  correlation  coefficient  for  the  failure  levels 
must  be  given  and  the  correlation  for  the  margin  of  survival  between  locations 
must  oe  determined. 

The  correlation  coefficient  for  the  margin  of  survival  at  any  two  loca¬ 
tions,  i  and  m,  can  be  computed: 


„  _  Cov  [Mt,  HJ 
i  ,m 

/VL'\J  V|MmJ 

Note  tnat  CovfM*,  rtj  -  E[Mj  -  E[Mt]  E[Mm] 

-  L,)(^  -  LJ]  -  E[MA]E[Mm] 


(77) 


now 


t^LXLmJ  *  p4,«/vtL*:]  v^Lm^  +  E^Lm^ 


wnere 


(>i  m  is  the  specified  correlation  coefficient  between  the  failure 
’  levels  at  locations  t  or  m 


t U  L  ]  *  £[Z  3  t[l  J 
l  m  l  m 


tC^zJ  3  £[L4]  E[Zm3  by  independence. 


Combining  the  above  gives: 


*  Z*(u)  Z^u)  +  £  £  Pjk  °k  ~ 

j*I  K=i  36k  (w }  36j  {u] 

♦  »tV,CLlJ  £CLm]  -  E[Zt]  £[LJ 


-  eczm]  eat]  -  ec«t]  tc-j 


wnich  contains  only  previously  determined  or  specified  elements. 

Failure  occurs  wnen  tne  margin  of  survival  is  negative,  and  the  probabil¬ 
ity  of  failure  can  therefore  oe  approximated  by 


Jf  3  p(,51  ^i  <  °)  s  Sp(Mi  <  0)  '  £D>0«i 

'  i=l  i*j 


<  0  Mj  <  0) 


+  ...  tp  Z  ",  <  o 


=P  I  ^  i 
\i*l 


wnere  in  is  the  number  of  failure  points  considered. 


«*.  •*. •*.  r_  . 
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This  series  may  be  approximated  by  truncating  the  series  of  summations 


after  several  terms.  The  problem  then  becomes 

f(us  Mi  <  o)  (80) 

where  s  is  a  subset  of  the  set  {1,  N}.  This  expression  can  be 

t  h 

evaluated  by  executing  the  n-fold  Integral  of  the  n  order  normal  pdf  of  Xj, 
n.  (Assume  that  the  set  s  contains  n  components  and  these  are 
denoted  X .,  j=l,  n.) 

J 

4  u 

The  n  order  normal  joint  pdf  is 

P„,,  ....  „n  (X,.  ....  X„)  . 

ei<p  (■  |  £  i;<xr“ t>  ru 

'  c  ial  jal  ' 

-  •  <  xi  <  •,  i*l,  n  (81) 

where  |S|  is  the  covariance  matrix  of  the  X.,  j«l,  n,  |S|  is  the 

k  U  k  U 

determinant  of  |S|,  r^j  is  the  cofactor  of  the  i  row,  j  column  element  in 

|S|. 

In  terms  of  Equation  81,  the  expression  80  is  written 


The  integral  can  be  modified  through  a  change  of  variables 


1*4 


1*1, 


to  obtain 


°i. 


[2*)nf2  | S| 1/2 


(83) 


(84) 


The  integral  can  be  further  modified  by  another  change  of  variables 


z,-  *  ,  i=l,  •••,  n 

yi 


(85) 


This  yields 


/  iz,  ...  /  iz„ 

\1£s  /  -<Vvi  (2 ^)n/2  |S| 


1/2  z\  •••  'z* 


n  n 


exp  ,  -  I  D  E  lillLll 
2  1-1  .1*1  Zi  Zj 


(86) 


I 


This  quantity  can  be  evaluatec  mmerically. 
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In  the  present  Investigation,  a  Simpson's  rule  scheme  Is  used  to  numeri¬ 
cally  evaluate  the  integral.  Each  of  the  integrals  evaluated  in  Equation  86 

is  an  element  in  one  of  the  sums  on  the  right  side  of  Equation  79. 

M 

Thus,  using  EfM^J,  VCM^],  and  Pt  m,  the  probability  of  failure  can  be 
obtained. 

d.  Numerical  example— A  nine  degree-of-freedoin  spring  mass  system  that 
models  a  structure  responding  only  axially  has  been  selected  to  illustrate  the 
analytical  procedure  (Fig.  18). 


Note:  All  masses  are  1000  kg. 


Figure  18.  Spring  mass  system. 


The  system  has  elasto-plastlc  springs.  Two  of  the  springs  near  the 

center  (springs  8  and  6)  exhibit  yield  at  their  mean  values  at  lower  levels 

tnan  do  tne  otners.  The  forcing  function  is  exerted  on  tne  first  mass  as 

-at 

shown  in  Figure  18  and  is  a  decaying  exponential  of  the  form  F  *  A  +  Be 
Three  of  the  parameters  that  descrioe  tne  system  and  two  that  describe  the 
applied  force  were  identified  as  random.  The  random  variaole  parameters  (mean 
and  variance)  and  tne  correlation  coefficients  between  sets  of  random 
variables  are  shown  in  Table  7. 
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The  problem  was  constructed  so  that  at  the  mean  values  of  the  parameters, 
the  maximum  strain  in  the  system  would  result  in  a  ductility  ratio  of  about 
10. 

The  partial  quotients  (that  estimate  tne  partial  derivatives  of  the 
response  measure  with  respect  to  each  of  the  random  parameters)  were  obtained 
by  a  series  of  deterministic  calculations  In  which  the  random  parameters  were 
varied.  The  deterministic  approach  was  the  numerical  solution  used  in  the 
previous  example. 


Random  Variables 


TABLE  7.  RANDOM  PARAMETERS  _ 


Mean  Variance  Correlation  Coefficients 


Spring  constant,, 
all  springs,  K 

1.0QE10 

1.00E18 

o 

• 

H 

c 

• 

M 

Yield  limit, 
springs  1-4,  7-9, 

1.00E08 

1.U0E14 

0.8 

Yield  limits, 
springs  5,  6,  Cj, 

8.O0E07 

7.23E13 

Peak  applied  force,  A+8 

1.00E08 

1.00E14 

Decay  coefficient,  a 

15 

2.25 

five  locations  were  selected  as  potential  failure  points;  one  near  the 
loading  (spring  2),  two  at  the  low-yield  springs  (springs  5  and  6),  and  two 
near  the  bottom  (springs  8  and  9).  The  failure  level  at  each  point  was 
considered  as  random;  the  means  and  the  variances  are  shown  in  Table  8.  The 
correlation  coefficients  Between  the  failure  levels  at  the  various  points  were 
also  estimated  and  are  shown  in  the  table. 


TntiLt  8.  FAILURE  LEVEL  MEAN  AND  VARIANCE 


Location 

Mean 

Variance 

Correlation 

Coefficients 

(Spring  dumber) 

5  6 

8 

9 

2 

0.06 

3.6E-05 

0.1  0.1 

0.3 

0.3 

5 

0.09 

8.16-05 

0.8 

0.1 

0.1 

6 

0.09 

8.1E-05 

0.1 

0.1 

8 

0.06  . 

3.5E-05 

0.7 

9 

0.06 

3.6E-05 

• 

Equations  74  and 

75  have  been 

useo  to 

estimate  the  mean 

and 

variance  of 

the  margin  of  survival  at  each  location.  Equation  78  estimates  the 
correlation  coefficients  between  the  margins  of  survival  at  the  various 
points.  This  information  has  been  used  as  input  in  an  approximation  of  the 
multivariate  normal  distribution,  and  a  numerical  form  of  Equation  79  has  been 
usea  to  compute  the  probaoility  of  system  failure.  The  analysis  shows  a 
failure  probability  of  56.4  percent. 


-v. : 


:_k 
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In  summary,  this  section  has  shown  the  development  of  an  approach  that 
predicts  the  reliability  of  MUF  systems  that  can  fail  at  any  of  a  number  of 
locations.  The  random  nature  of  the  structural  system  and  the  load  have  been 
addressed.  In  addition,  the  failure  level  has  been  considered  as  a  random 
property.  The  development  is  general  in  that  the  random  parameters  are 
treated  as  correlated  rather  than  Independent. 


IV.  DEVELOPMENT  OF  A  STOCHASTIC  FINITE  ELEMENT  COOE 


1.  INTRODUCTION 

The  previous  section  addressed  the  use  of  a  deterministic  finite  element 
code  to  calculate  tne  mean  and  to  approximate  the  variance  of  selected 
response  measures  for  nonlinear  problems  with  random  input  and  material  pro¬ 
perties.  The  approach  required  repetitive  calculations  to  estimate  tne  mean 
ana  tne  behavior  near  the  mean  so  that  the  variance  could  be  estimated.  Since 
nonlinear,  dynamic  problems  are  random  processes  it  would  be  more  useful  to 
nave  an  approacn  tnat  tracxs  tne  evolution  of  tne  random  behavior  of  response 
measures  witn  time.  This  section  addresses  such  an  approacn.  The  equations 
in  a  typical  numerical  approach  may  be  used  to  compute  the  variance  and  covar¬ 
iance  of  various  response  parameters;  but  this  section  carries  the  results 
only  as  far  as  the  random  description  of  the  response  measures.  The  extension 
for  proDabilistic  failure  analysis  is  the  same  as  for  the  previous  section. 

2.  THEORETICAL  DEVELOPMENT 

Consiaer  a  one-dimensional  nonlinear  finite  element  algoritnm  which 
proceeds  as  follows.  Tne  acceleration  at  location  i  and  at  tne  jtn  time  step 
is 
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where 


a  *  acceleration 
F  »  applied  force 
FI  »  internal  force 
m  *  mass 

superscript  i  is  a  spatial  index 
subscript  j  is  a  time  index 

The  velocity,  V,  at  location  i,  time  step  j,  for  a  At  time  step  is 

V  *  V1  +  a ’at 
J  J-1  J 

and  the  displacement,  d,  becomes 

d^  *  d*  +  At 
J  J-1  j 

Using  a  simple  definition  of  strain  based  on  the  original  length,  i,  the 
strain  e  becomes 


and  stress,  o,  is  addressed  in  some  nonlinear,  inelastic  relationship  to 
strain 


(v 


S-l’  0n,j-l’ 


•) 
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where  em  and  0mj_j  are  memory  parameters  (from  the  last  time  step)  that  are 
usually  required  for  path  dependent  models.  E  and  other  parameters  are  the 
moduli  and  strains  that  are  used  In  the  functional  form  of  the  nonliner 
model . 

An  example  of  a  nonlinear  model  (used  in  subsequent  calculations  in  this 
report)  is  a  two-modulus  elastic-plastic  model  that  unloads  on  the  elastic 
modulus.  For  strains  less  than  some  specified  strain  ej,  the  material  model 
is  elastic  with  modulus  t,  that  is 


wnen  strains  lie  in  the  range  <  e\  <  e2,  then  a  second  modulus,  £2,  becomes 

J 

effective  and 

^  •  &i  *  (ej  -  ei)E2 

and  when  tne  strain  exceeds  tl  the  material  is  plastic  with 
a’  *  cei  +  (el-e2)E2 

Tii i s  particular  model  can  treat  unloading  differently  than  loading  and  thus 
dissipate  energy.  This  treatment  requires  memory  parameters.  The  memory 
parameters  are  the  maximun  stress  and  strain  attained,  so  that  at  unload  the 
stress  can  be  computed  from 


The  last  computation  for  a  time  step  in  this  finite  element  process  is 
the  calculation  of  the  internal  force  FI  for  the  next  time  step. 

fi3  ■  W"1  -  t,J)A 

where 

A  »  cross  section  area  of  the  element 

o°.  »  0 
J 


After  all  elements,  i  *  1.  .  .N,  and  nodes,  i  *  1 — N  +  1,  are  updated,  the 
information  generated  at  j  is  used  to  compute  the  response  at  j  +  1. 


The  fact  tnac  all  the  memory  required  for  calculation  of  the  response 
measures  for  the  current  time  step  lies  in  the  previous  time  step  information 
indicates  that  this  may  be  a  Markov  Process.  A  Markov  Process  is  one  in  whicn 
the  properties  of  a  random  process  can  be  computed  given  the  state  at  a  pre¬ 
vious  time  and  the  time  since  that  state.  The  random  properties  of  each  of 
the  response  measure  parameters  can  be  computed  from  the  previous  time  step 
data.  For  example  the  first  output  measure  of  interest  that  is  computed  by 


the  deterministic  algorithm  is  the  velocity 

vj  *  Vj-1  +  ajAt 

but 


so  that 


so  that 


Var 


M]  ■  var[vj.j  - 2  ^  v*rtFj]  * var  [Fri-i] 

m  m 


It  Is  apparent  that  It  will  be  necessary  to  compute 


Var 


[v]FIj] 


,1  • 


for  the  next  time  step  and  that  will  be  calculated  when  FI j  is  found. 
The  displaced  mean  and  variance  are 


*  deterministic  calculations  using  mean  values 


-K]  ■  fK2]  • 

E[d]2]  -  V.r[d].J  *  2«Co» 

*  *«{*]] 


2 


Cov  d 


i 

j-1 


The  equations  for  the  random  properties  for  the  remaining  parameters  In  the 
one-dimensional  finite  element  algorithm  have  been  developed  and  are  available 
in  the  computer  program  found  in  the  Appendix.  The  program  computes  the  evolu 
tlon  of  the  multivariate  random  distribution  of  the  response  measures,  includ¬ 
ing  mean  and  variance  for  each  parameter  and  covariance  between  parameters. 


[•] 


how,  if  a  norma)  or  a  lognormal  distribution  is  selected  to  represent  the 


form  of  tne  distribution,  then  the  parameters  exist  to  completely  define  the 
distribution  at  any  time  step. 

fcven  tnougn  tne  response  measures  form  a  multivariate  distribution  and 
integration  would  require  an  n-fold  integration  scheme,  the  marginal  prob¬ 
ability  density  functions  can  be  assumed  to  be  normal  or  lognormal  if  the 
joint  probability  density  is  likewise  normal  or  lognormal,  so  that  some 
information  can  be  obtained  from  the  marginal  density  functions. 

3.  NUrttKlCAL  tXAMPLt 

A  bar  5  m  long,  with  properties  as  shown  in  Table  9,  is  subjected  to  a 
decaying  exponential  load  of  the  form 


where 

m  »  4.15  MPa 

a  *  1000 

t  *  time  in  seconds 

The  above  problem  has  been  solved  using  a  finite  element  procedure  witn  A 
as  a  ranaom  parameter  with  mean  4.15  and  standard  deviation  0.4.  Response 
measures  of  the  first  two  nodes  are  shown  in  Figures  19  through  22.  A  one- 
standard  deviation  bound  on  the  mean  is  shown  for  each  parameter.  A  second 
problem  is  solved  whicn  included  random  material  properties  as  shown  in  Table 
10.  Tnis  solution  is  also  shown  on  Figures  19  through  22. 


TABLE  9.  PROPERTIES  OF  BAR 


Area  *  O.Ol  m2 
Density  =  7840  kg/m3 
*  2  E+ll  Pa 
E2  *  1  E+ll  Pa 


«1 

*2 


-0.022 

-0.005 


where  E1 ,  E2,  e  ,  and  e2  are  defined 
by  the  following: 


Strain  (compression) 

The  model  unloads  on  the  initial  modules,  E 


TABLE  10. 

RANDOM  MATERIAL 

PROPERTIES 

Standard 

Parameter 

Mean,  Pa 

deviation,  1 

El 

2  E+ll 

2  E+10 

1  E+ll 

1  E+10 

el 

-0.002 

0.0002 

€2 

-0.005 

0.0001 
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tjCmJL 


Velocity,  m/s  .  ,  .  Velocity,  »/s 


rain 


Time,  s 


Stress.  MPa  Stress,  MPa 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


1.  CONCLUSIONS 

Procedures  have  been  developed  to  compute  the  probability  of  failure  of 
structures  subjected  to  blast  and  shock  loads.  The  requirements  are  as 
follows: 


a.  The  characteristics  of  the  random  load  must  be  specified. 

b.  The  structure  system  parameters  tnat  are  considered  random  must  be 
adequately  specified. 

c.  ueterministic  tools  to  compute  the  response  must  be  availaDle. 

a.  The  random  nature  of  the  failure  criteria  must  be  characterized. 

There  are  limitations  to  the  presented  approaches.  For  the  approach 
presented  for  S(/F  systems,  the  computational  storage  requirement  is  quite 
large.  Methods  to  reduce  the  requirements  rely  on  narrow  band  response 
assumptions  wnicn  are  not  correct  for  hlgnly  nonlinear  systems.  The 

f 

discretization  of  response  space  Introduces  further  numerical  error  and  must 
oe  carefully  selected  with  respect  to  time  step  and  response  measure  nistory 
results. 

The  one  approach  presented  for  MDF  representations  is  limited  by  the 
requirement  to  truncate  the  Taylor  series  at  a  reasonable  order  so  as  to  limit 
tne  number  of  deterministic  calculations  that  mu^t  be  run  to  provide  the 
numerical  approximations  of  the  partial  differences.  This  truncation  has  been 
at  the  linear  terms  for  the  example  problems,  and  is  adequate  as  long  as  the 


distribution  of  the  response  about  Its  mean  Is  narrow  enough  so  that  the 
linear  approximation  Is  adequate.  Otherwise  higher  order  terms  must  be  main¬ 
tained.  Tne  linear  terms  can  be  adequate  for  nonlinear  response  If  the  peak 
range  over  which  the  linear  approximation  Is  used  accounts  for  a  significant 
portion  of  tne  peak  response  distribution.  That  Is,  for  reasonable  bounds  on 
the  peax  response  distribution,  the  peak  response  function  must  be  linear  in 
the  random  parameters. 

Tne  second  NUF  approach  (finite  element)  would  require  extensive  develop¬ 
ment  for  two-dimensional  problems. 

The  approacnes  do  not  address  model  uncertainty.  Model  uncertainty 
snould  not  be  addressed  witn  a  random  property  approach  and  appropriate  pro¬ 
cedures  are  outside  the  scope  of  this  report. 

In  practical  application  of  the  approaches  developed  in  this  investiga¬ 
tion,  the  most  severe  limitation  Is  an  inadequate  data  base  from  which  to 
establisn  random  character  of  the  loads,  the  structure  system  parameters,  and 
the  failure  criteria.  Testing  is  very  expensive  because  of  the  nature  of  the 
loads;  ana  reasonable  tests  provide  some  damage  to  the  structure,  making  them 
difficult  to  repeat.  Further  definition  of  what  can  be  correctly  addressed  as 
a  random  variaole  is  needed. 

2.  EXTENSIONS 

Tne  most  useful  steps  beyond  this  study  include: 

a.  more  thorough  criteria  for  determining  the  order  of  truncation  of  the 


Taylor  series.  This  would  be  developed  mostly  by  extensive  experience  with  a 
variety  of  problems. 


b.  The  one-dimensional  finite  element  model  could  be  expanded  to  two- 
dimension.  It  may  be  necessary  to  approximate  some  of  the  relationships  to 
keep  the  development  to  a  manageable  state.  More  experience  with  the  one- 
dimensional  model  would  be  helpful. 

c.  Tne  distribution  of  response  measures  could  be  used  to  validate 
whetner  a  deterministic  model  could  produce  random  response  for  random  Inputs 
In  the  same  way  as  it  does  for  the  real  system.  For  example,  given  a  reason¬ 
able  random  Input  loading,  both  a  mathematical  model  and  a  test  could  be 
conducted  for  concrete  specimens.  The  distribution  of  the  strain  from  both 
could  De  compared  to  validate  tne  mathematical  model.  Because  of  tne  complex 
stress/strain  paths  that  are  investigated,  normalization  of  data  along  path 
segments  would  probably  be  required. 
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